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ABSTRACT 

This  investigation  is  directed  toward  understanding  the  role  of  coupled  me¬ 
chanical  and  thermal  effects  in  the  linear  stability  of  an  isothermal  antiplane  shear 
motion  which  involves  a  single  planar  phase  boundary  in  a  non- elliptic  thermoe¬ 
lastic  material  which  has  multiple  elliptic  phases.  When  the  relevant  process  is 
static — so  that  the  phase  boundary  does  not  move  prior  to  the  imposition  of  the 
disturbance — it  is  shown  to  be  linearly  stable.  However,  when  the  process  involves 
a  steadily  propagating  phase  boundary  it  may  be  linearly  unstable.  Various  con¬ 
ditions  sufficient  to  guarantee  the  linear  instability  of  the  process  are  obtained. 
These  conditions  depend  on  the  monotonicity  of  the  kinetic  response  function — a 
constitutively  supplied  entity  which  relates  the  driving  traction  acting  on  a  phase 
boundary  to  the  local  absolute  temperature  and  the  normal  velocity  of  the  phase 
boundary — and,  in  certain  cases,  on  the  spectrum  of  wave- numbers  associated 
with  the  perturbation  to  which  the  process  is  subjected.  Inertia  is  found  to  play 
an  insignificant  role  in  the  qualitative  features  of  the  aforementioned  sufficient 
conditions.  It  is  shown,  in  particular,  that  instability  can  arise  even  when  the 
normal  velocity  of  the  phase  boundary  is  an  increasing  function  of  the  driving 
traction  if  the  temperature  dependence  in  the  kinetic  response  function  is  of  a 
suitable  nature.  Significantly,  the  instability  which  is  present  in  this  setting  oc¬ 
curs  only  in  the  long  waves  of  the  Fourier  decomposition  of  the  moving  phase 
boundary,  implying  that  the  interface  prefers  to  be  highly  wrinkled. 
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1.  INTRODUCTION 

Recently  [9],  motivated  by  a  desire  to  determine  whether  continuum  me¬ 
chanical  models  for  displacive  solid-solid  phase  transformations  can  predict  the 
emergence  of  plate-like  or  dendritic  structures  from  states  involving  planar  phase 
boundaries,  a  purely  mechanical  two- phase  dynamical  process  in  a  non- elliptic 
generalized  neo-Hookean  material  was  considered.  The  process  involved  an  an- 
tiplane  shear  motion  with  a  single  steadily  propagating  planar  phase  boundary 
separating  high  and  low  strain  elliptic  phases  of  the  relevant  material.  In  a  frame 
moving  with  the  phase  boundary,  the  shear  strain  field  was  piecewise  homogeneous 
and  the  angle  between  the  limiting  values  of  the  gradient  of  the  out-of-plane  dis¬ 
placement  field  on  either  side  of  the  phase  boundary  was  zero — so  that  the  phase 
boundary  was,  for  each  instant  of  the  motion,  of  normal  type.  The  linear  stability 
of  this  process  with  respect  to  a  broad  class  of  perturbations  was  then  investi¬ 
gated.  It  was  shown  that  a  necessary  and  sufficient  condition  for  the  process  to 
be  linearly  stable  was  that  the  kinetic  response  function — which  gives  the  driving 
traction  acting  on  a  phase  boundary  in  terms  of  the  normal  velocity  of  the  phase 
boundary,  or  vice-versa — be  a  locally  increasing  function  of  its  argument  at  the 
value  corresponding  to  the  base  process.  A  necessary  consequence  of  this  stability 
criterion  is  that,  in  order  for  the  process  to  be  unstable,  the  kinetic  response  func¬ 
tion  must  exhibit  a  non-monotonic  dependence  on  its  argument.  Non-monotonic 
kinetic  response  functions  are  admissible  under  the  Clausius- Duhem  version  of 
the  second  law  of  thermodynamics  (specialized  to  isothermal  conditions  for  the 
purposes  of  the  purely  mechanical  process  discussed  in  [9]);  the  work  of  Owen, 
Schoen  &  Srinivasan  [15]  implies,  furthermore,  that  a  non-monotonic  relation 
between  interfacial  driving  traction  and  normal  velocity  may  be  operative  in  the 
unstable  kinetics  which  are  observed  to  accompany  the  emergence  and  growth  of 
plate-like  structures.  Under  such  kinetics,  the  results  obtained  [9]  suggest  than 
an  evolution  from  a  planar  to  a  plate-like  phase  boundary  morphology  might  be 
possible  with  the  confines  of  a  purely  mechanical  theory. 
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Thermal  effects  are  manifestly  absent  from  the  purely  mechanical  investi¬ 
gation  in  [9].  The  experimental  work  of  Clapp  &  Yu  [5],  Grujicic,  Olson  & 
Owen  [10]  and  Cong  Dahn,  Morphy  &  Rajan  [6]  indicates  that  temperature 
effects  do  play  an  intrinsic,  if  not  entirely  understood,  part  in  the  kinetics  of 
phase  boundaries  in  displacive  solid-solid  phase  transformations.  The  investiga¬ 
tion  which  follows  is,  therefore,  directed  toward  understanding  the  outcome,  with 
regard  to  the  morphological  stability  of  states  involving  planar  phase  boundaries, 
when  thermal  as  well  as  mechanical  effects  are  taken  into  consideration  in  a  model 
for  displacive  solid-solid  phase  transformations.  Of  particular  interest  is  the  ques¬ 
tion  of  whether  thermal  effects  allow  for  an  evolution  from  planar  to  plate-like 
phase  boundary  morphology  under  kinetics  which  are  mechanically  stable  in  the 
sense  of  [9].  The  paper  is  organized  as  follows. 

Chapter  2  is  dedicated  to  preliminaries.  Following  a  synopsis  of  the  notation 
to  be  used,  Section  2.1  introduces  the  kinematics  and  fundamental  balance  prin¬ 
ciples  which  will  be  needed  thereafter.  Section  2.2  focuses  on  the  rate  of  entropy 
production  due  to  the  presence  of  phase  boundaries  and  introduces  the  associ¬ 
ated  notion  of  the  driving  traction  acting  on  a  phase  boundary.  In  Section  2.3  a 
thermoelastic  material  is  defined  and  in  Section  2.4  the  particular  class  of  ther¬ 
moelastic  materials  which  will  be  used  in  the  forthcoming  analysis  is  introduced. 
Section  2.5  is  concerned  with  the  kinetic  relation  and  allied  kinetic  response  func¬ 
tion.  In  the  final  section  of  Chapter  2  the  kinematics  are  specialized  to  those  of 
antiplane  shear  and  a  thermoelastic  antiplane  shear  motion  is  defined. 

Chapter  3  is  devoted  to  the  linear  stability  analysis  of  an  isothermal  two- 
phase  process  which  involves  a  steadily  propagating  planar  phase  boundary  in 
an  arbitrary  thermoelastic  material  within  the  class  introduced  in  Section  3.4. 
The  relevant  process,  which  is  a  straightforward  generalization  of  that  used  in  the 
purely  mechanical  investigation  [9],  is  introduced  in  Section  3.1,  while  the  class 
of  perturbations  to  which  it  will  be  subjected  is  put  forth  in  Section  3.2.  Each 
admissible  perturbation  involves,  in  general,  a  disturbance  of  the  configuration 


of  the  phase  boundary  and  of  the  displacement,  velocity  and  temperature  fields 
in  a  small  neighborhood  of  the  phase  boundary.  All  disturbances  are  assumed 
to  be  small  in  some  appropriate  sense.  The  kinematics  of  the  perturbation  are 
restricted  to  those  of  antiplane  shear.  It  is  assumed  that  the  post-perturbation 
process  is  a  thermoelastic  antiplane  shear  and  involves  only  one  phase  boundary. 
Sections  3.3  and  3.4  address,  respectively,  the  linearization — about  the  unper¬ 
turbed  process — of  the  field  equations  which  hold  away  from  the  phase  boundary 
and  the  jump  conditions  and  kinetic  relation  which  hold  on  the  phase  boundary. 
After  a  specialization  of  the  base  process,  a  summary  of  the  complete  linearized 
system  of  field  equations,  jump  conditions,  kinetic  relation  and  boundary  and 
initial  conditions  which  describe  the  process  generated  by  the  perturbation  is 
presented  in  Section  3.5.  As  in  [9],  both  the  inertial  and  inertia-free  cases  are 
included.  The  combined  results  of  Sections  3.6  and  3.7  show  that  whenever  it 
is  static,  regardless  of  the  presence  of  inertial  effects,  the  base  process  is  linearly 
stable  with  respect  to  all  perturbations  of  the  type  introduced  in  Section  3.2. 
Section  3.8  deals  with  the  case  where  the  base  process  involves  an  interface  prop¬ 
agating  at  non-zero  velocity.  A  normal  mode  analysis  is  performed  which  leads 
to  a  variety  of  conditions  sufficient  for  the  instability  of  the  undisturbed  process. 
These  conditions  depend  on  the  monotonicity  properties  of  the  kinetic  response 
function.  Highlighted  in  Section  3.9  is  one  set  of  sufficient  conditions  which  is 
of  particular  interest.  The  relevant  conditions  alter  the  conclusions  reached  in 
the  purely  mechanical  context  considered  in  [9]  in  two  ways.  First,  in  contrast 
to  the  results  obtained  in  the  latter  setting,  instability  may  arise  even  when  the 
normal  velocity  of  the  phase  boundary  is  a  monotonically  increasing  function  of 
driving  traction  as  long  as  the  temperature  dependence  in  the  kinetic  response 
function  is  of  an  appropriate  nature.  Second,  the  instability  that  arises  in  these 
thermomechanical  circumstances  occurs  only  in  the  long  waves  of  the  the  Fourier 
decomposition  of  the  moving  phase  boundary,  suggesting  that  the  interface  favors 
a  highly  wrinkled  configuration.  This  conclusion  is  akin  to  that  reached  in  simi- 
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lar  linear  stability  analyses  performed  within  the  context  of  models  for  dendritic 
crystal  growth  where  an  otherwise  thermally  unstable  process  is  stabilized  for 
sufficiently  large  wave-numbers  by  the  inclusion  of  surface  tension  at  the  inter¬ 
face  in  lieu  of  including  mechanical  effects.1  The  final  topic  addressed  in  Section 
3.9  pertains  to  the  physical  suitability  of  kinetic  response  functions  which  are 
mechanically  stable  but  thermally  unstable. 


1  See,  for  example,  Langer  [13],  Mullins  &  Sekerka  [14]  and  Strain  [16]. 
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2.  PRELIMINARIES 

2.1.  Notation,  kinematics  and  balance  principles.  In  the  following  ]R 
and  (F  denote  the  sets  of  real  and  complex  numbers.  The  intervals  (0,  oo)  and 
[0,  oo)  are  represented  by  jR+  and  M+.  The  symbol  2R",  with  n  equal  to  2  or  3, 
represents  real  n-dimensional  space  equipped  with  the  standard  Euclidean  norm. 

—  O 

If  U  is  a  set,  then  its  closure,  interior  and  boundary  are  designated  by  U,  U ,  and 
dU,  respectively.  The  complement  of  a  set  V  with  respect  to  U  is  written  as 
U  \  V.  Given  a  function  ip  :  U  — +  W  and  a  subset  V  of  U,  ip(V)  stands  for  the 
image  of  V  under  the  map  ip. 

Vectors  and  linear  transformations  from  JR3  to  JR3  (referred  to  herein  as 

tensors )  are  distinguished  from  scalars  with  the  aid  of  boldface  type — lower  and 

upper  case  for  vectors  and  tensors,  respectively.  Let  a  and  b  be  vectors  in  M3, 

their  inner  product  is  then  written  as  a  •  b;  the  Euclidean  norm  of  a  is,  further, 

written  as  |a|  =  >/a  •  a.  The  set  of  unit  vectors — that  is,  vectors  with  unit 

Euclidean  norm — in  M3  is  designated  by  Af.  The  symbol  C  refers  to  the  set 

+ 

of  tensors,  C+  denotes  the  set  of  all  tensors  with  positive  determinant,  and  S 
stands  for  the  collection  of  all  symmetric  positive  definite  tensors.  If  F  is  in  C 
then  Fr  represents  its  transpose;  if,  moreover,  det  F  /  0,  then  the  inverse  of  F 
and  its  transpose  are  written  as  F-1  and  F~r,  respectively.  The  notation  a®  b 
refers  to  the  tensor  A,  formed  by  the  outer  product  of  a  with  b,  defined  such  that 
Ac  =  (b  •  c)a  for  any  vector  c  in  1R3.  If  A  and  B  are  tensors  then  their  inner 
product  is  written  as  A  •  B  =  tr  ABT. 

When  component  notation  is  used,  Greek  indices  range  only  over  {1,2}; 
summation  of  repeated  indices  over  the  appropriate  range  is  implicit.  A  subscript 
preceded  by  a  comma  denotes  partial  differentiation  with  respect  to  the  corre¬ 
sponding  coordinate.  Also,  a  superposed  dot  signifies  partied  differentiation  with 
respect  to  time. 

Consider  now  a  body  B  which,  in  a  reference  configuration,  occupies  a  region 
%  contained  in  1R3.  A  motion  of  B  on  a  time  interval  T  C  1R  is  cheiracterized  by 
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a  one-parameter  family  of  invertible  mappings  y(-,  t)  :  TZ  — *  Ht,  with 

y(x,t)  =  x  +  u(x,<)  V(x,i)eM,  (2.1.1) 

where  M  —  Tt  x  T  represents  the  trajectory  of  the  motion.  Assume  that  the 
deformation  y,  or  equivalently  the  displacement  u,  is  continuous  and  possesses 
piecewise  continuous  first  and  second  partial  derivatives  on  M.  Let  St  be  the  set 
of  points  contained  in  1Z  defined  so  that,  at  each  instant  t  in  T,  y(-,t)  is  twice 
continuously  differentiable  on  the  set  7 Z\  St-  Let  the  set  E  be  defined  by 

E  =  {(x,t)|x  6  St,t  €  T).  (2.1.2) 

Introduce  the  deformation  gradient  tensor  F  :  M  \  E  — »  C  by 

F(x,<)  =  Vy(x,t)  V(x,t)  £  M\E,  (2.1.3) 

where  the  associated  Jacobian  determinant,  J  :  M  \  E  — ►  1R,  of  y  is  restricted  to 
be  strictly  positive  on  its  domain  of  definition: 

J(x,<)  =  det F(x,f)  >  0  V(x,t)eM\E. 

+ 

Hence,  F  :  M  \  E  — ►  £+.  The  left  Cauchy-Green  tensor  G  :  M\E  —*  S 
corresponding  to  the  deformation  y  is  given  by 

G(x, t)  =  F(x, t) Ft(x, t)  V(x, t)  €  M\E.  (2.1.4) 

The  deformation  invariants  associated  with  y  exist  on  M  \  E  and  are  supplied 
through  the  fundamental  scalar  invariants  of  G: 


J,(G)=trG,  J2(G)  =  i((trG)2-tr(G2)),  J3(G)  =  detG.  (2.1.5) 
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With  the  above  kinematic  antecedents  in  place  introduce  the  nominal  mass 
density  p  :  72.  — ►  ]R+,  the  nominal  body  force  per  unit  mass  b  :  M  — »  IF?3,  and  the 
nominal  stress  tensor  S  :  M.  \  £  — ►  C,  and  suppose  that  p  is  constant  on  72  and  b 
is  continuous  on  M.,  while  S  is  piecewise  continuous  on  M,  continuous  on  M\  H, 
and  has  a  piecewise  continuous  gradient  on  M.  Let  p„  be  the  mass  density  in 
the  deformed  configuration  associated  with  y.  Given  a  regular  subregion  V  of 
72,  with  dV  n  St  a  set  of  measure  zero  in  dV  for  each  t  in  T,  let  m  :  dV  — ►  Af 
denote  the  unit  outward  normal  to  dV.  Then  the  global  balance  laws  of  mass, 
linear  momentum,  and  angular  momentum  require  that 


on  T, 


(2.1.6) 


J  Sm<L4  +  J  pbdV  =  J  pixdV  on  T, 

dV  V  V 


(2.1.7) 


and 


JyASmdA  +  JyApbdV  =  jyf\p\xdV  on  T,  (2.1.8) 
dV  v  v 


respectively,  for  every  such  regular  subregion  V  contained  in  72.. 

Next,  introduce  the  nominal  internal  energy  per  unit  mass  e  :  M  \  £  — ►  JR, 
the  nominal  heat  flux  q  :  M  \  S  — *  M3,  and  the  nominal  heat  supply  per  unit 
mass  r  :  M\S  — ►  1R.  Suppose  that  e  and  q  are  piecewise  continuous  on  M , 
continuous  on  M  \  E,  and  have  piecewise  continuous  partial  derivatives  on  M, 
and  that  r  is  continuous  on  M .  The  first  law  of  thermodynamics  requires  that 


•  u  +  qm 


v 


p(b-  u  +  r)dV 


J  p{e  +  ±\u\*)dV  on  T,  (2.1.9) 
v 


for  every  regular  subregion  V  contained  in  72  such  that  dV  PI  St  a  set  of  measure 
zero  in  dV  for  each  t  in  T. 


Finally,  introduce  the  nominal  entropy  per  unit  mass  r)  :  M\E  —*  JR  and 
the  nominal  absolute  temperature  6  :  M.  —*  JR.  Stipulate  that  77  is  piecewise 
continuous  on  M,  continuous  on  M  \  E,  and  has  piecewise  continuous  first  partial 
derivatives  on  M,  and  that  9  is  continuous  on  M  with  piecewise  continuous  first 
partial  derivatives  on  M.  The  Clausius- Duhem  version  of  the  second  law  of 
thermodynamics  requires  that  the  rate  of  entropy  production  F(-;  V)  :  T  — ►  JR 
satisfies 

/’(•;*>)  =  J  pr,dv- J^dA- J  ^dV>0  on  T,  (2.  MO) 

V  dV  V 


for  every  regular  subregion  V  contained  in  R  such  that  dV  fl  St  a  set  of  measure 
zero  in  dV  for  each  t  in  T. 

Localization  of  the  balance  laws  (2.1.6)-(2.1.9)  and  the  imbalance  law  (2.1.10) 
at  an  arbitrary  point  contained  in  the  interior  of  M  \  E  yields  the  following 
familiar  field  equations  and  field  inequality: 


p = P.(?)J 

on 

■M\r, 

V  •  S  +  pb  =  pu 

on 

M\E, 

SFr=  FSr 

on 

M\E, 

S-F  +  V-  q  +  pr  =  pe 

on 

M\E. 

on 

M\E. 

(2.1.11) 


Suppose,  from  now  on,  that  the  set  St  is  a  regular  surface  for  every  t  in  T. 
The  set  E  then  represents  the  trajectory  of  a  surface  of  discontinuity  in  F,  S  and, 
perhaps,  e,  q  and  tj.  Let  g(-,t)  denote  a  generic  field  quantity  g(-,t)  :  St  —*  JR 
which  is  discontinuous  across  St  at  the  instant  t  in  T.  Define  the  jump  [<?(•,<)] 
of  g(-,t)  across  St  by 


where  n(-,t)  :  S<  — *  is  a  unit  normal  to  St  at  each  t  in  T.  Then,  localization  of 
(2.1.6)-(2.1.10)  at  an  arbitrary  element  of  E  yields  the  following  jump  conditions 


[P*(y)^l  =  0  on  E, 

[SnJ  +  pVnful  =  0  on  E, 

,  ,  (2.1.13) 

[Sn  •  u]  +  +  ||u|2]  +  [q  ■  n]  =  0  on  E, 

pVnlnl  +  ^[q  •  n]  <  o  on  E, 

where  Vn(^t)  :  St  — *  1R  is  the  component  of  the  velocity  of  the  surface  St  in  the 
direction  of  n(-,<)  at  the  instant  t  in  T. 

Equations  (2.1.11)i  and  (2.1.13)i  are,  evidently,  completely  decoupled  from 
equations  (2.1.11)2,3,4,5  and  (2.1.13)2,3,45  that  is,  given  a  solution  to,  say,  a  bound¬ 
ary  value  problem  involving  (2.1.11)2,3,4,5  and  (2.1.13)2,3,4,  p *  can  be  calculated 
a  posteriori.  For  this  reason  equations  (2.1.11)i  and  (2.1.13)i  will  be  disregarded 
in  the  subsequent  analysis. 

In  this  investigation  an  inertia-free  motion  is  defined  as  one  wherein  the 
inertial  terms  on  the  right  hand  sides  of  the  global  balance  equations  (2.1.7)  and 
(2.1.8)  are  replaced  by  the  zero  vector.  In  the  context  of  an  inertia-free  motion 
the  field  equation  (2.1.11)2  simplifies  to  read 

V  •  S  +  pb  =  0  on  M\E ,  (2.1.14) 

and  the  jump  condition  (2.1.13)2  becomes 

fSnJ  =  0  on  E.  (2.1.15) 

Equations  (2.1.11)1,3,4,5  and  (2.1.13)i,3,4  remain,  of  course,  unaltered. 

In  addition  to  the  jump  conditions  given  in  (2.1.13)  in  the  inertial  case  or 
(2.1.13)i,3,4  and  (2.1.15)  in  the  inertia-free  case,  the  stipulated  continuity  of  y 


and  6  gives  the  following  kinematic  jump  conditions 


[u|  =  0  on  E,  |[0]  =  O  on  E.  (2.1.16) 

2.2.  Rate  of  entropy  production  and  driving  traction.  Using  the  field 
equations  (2.1.11),  the  jump  conditions  (2.1.13)  ,  and  the  assumed  smoothness  of 
the  deformation  y,  Abeyaratne  &  Knowles  [1]  have  demonstrated  that  for  any 
continuum  the  rate  of  entropy  production  /'(•;  V)  can,  for  any  regular  region  V 
contained  in  71,  be  represented  in  the  form 

m  V)  =  r,oc(t;  V)  +  rcon(t;  V)  +  r,(t;  V)  Vt  G  T,  (2.2.1) 

where  iloc(-;"P),  rcon(-;  V),  and  ra(-;  V)  are  defined  by 

flo.(-;P)=  j  i(S  -  F  -  +  «7,))dV  on  T, 

v\st 

rco„(  ;P)=  /  veiv  on  T,  (2.2.2) 

v\s, 

r.(-;P)=  J  i(pM-«S)).[F])KdV  on  T, 

Vn  S, 

with  rp  :  M  \  E  — *  JR  representing  the  nominal  Helmholtz  free  energy  per  unit 
mass  in  defined  in  terms  of  e,  9  and  rj  by 

xp  =  e  —  Brj  on  M\E,  (2.2.3) 

and,  where — given  a  generic  field  quantity  g(-,t)  :  St  —*  JR  which  jumps  across  St 
at  the  instant  t  in  T — ((</(-,<)))  is  defined  through 
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The  representation  (2.2.1)  additively  decomposes  the  total  rate  of  entropy 
production  r(-;'P),  at  the  instant  t  in  T,  in  the  regular  region  V  contained  in  % 
into  three  parts.  The  first  two  terms  in  the  decomposition  floc(s  V)  and  rco „(•;  V) 
are  the  contributions  to  the  rate  of  entropy  production  due,  respectively,  to  local 
mechanical  dissipation  and  heat  conduction  away  from  the  surface  St,  while  the 
third  term  ra(-\V)  represents  the  entropy  production  rate  due  to  the  motion  of 
the  surface  St. 

Motivated  by  (2.2.2)3  define  the  driving  traction  /(•,<)  :  5<  — ►  1R  which  acts 
on  the  surface  St  at  the  instant  t  in  T  by 

/(•■<)  «/W-,*)]-«S(-,0»-lF(-,01  on  5‘  V,6T-  <2-2-5) 

In  the  absence  of  inertial  effects  it  can  be  demonstrated  that  (2.2.5)  reduces  to 

/(•.«)  =  on  5,  VfeT,  (2.2.6) 

+  — 

where  S(-,t)  (reap.,  S (•,<))  is  the  limiting  value  of  the  field  S(-,t)  on  the  side  of 
the  interface  into  which  the  unit  normal  n (•,<)  is  (reap.,  is  not)  directed  at  the 
instant  t  in  T. 

Now,  from  (2.2.1)  and  (2.2. 2)3,  localization  of  the  imbalance  law  (2.1.10)  at 
an  arbitrary  element  of  E  yields  the  following  alternative  to  (2.1.13)4: 

fVn>  0  on  E,  (2.2.7) 

with  /  given  by  (2.2.5)  or  (2.2.6)  depending  on  whether  inertial  effects  are  included 
or  not.  Observe  ,  from  (2.2.2),  that  under  isothermal  conditions  the  total  rate  of 
entropy  production  /’(•;  V)  in  a  region  V  takes  the  form  of  the  rate  of  mechanical 
dissipation  per  unit  temperature. 

2.3.  Finite  thermoelasticity.  Let  B  be  composed  of  a  homgeneous 
thermoelastic  material.  Then  there  exists  a  Helmholtz  free  energy  potential 
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rp  :  £+  x  2R+  — ►  2R  such  that  the  nominal  Helmholtz  free  energy  per  unit  mass 
rp,  the  nominal  stress  tensor  S,  and  the  nominal  entropy  per  unit  mass  are  given 
in  terms  of  rp  as  follows: 

xp  =  xp(  F,0)  on  M\E, 

S  =  (4f(F,9)  on  M\E,  (2.3.1) 

tj  =  —xpo(F,  0)  on  M  \  E. 

It  is  assumed  that  xp  is  once  continuously  differentiable  and  piecewise  twice  contin¬ 
uously  differentiable  on  £+  x  1R+.  The  nominal  heat  flux  q  is,  for  a  thermoelastic 
material,  given  by  a  heat  flux  response  function  q  :  C+  x  JR+  x  Ft3  — ►  1R 3  so  that 

q  =  q(F,0,V0)  on  M\E.  (2.3.2) 

It  is  assumed  that  q  is  piecewise  twice  continuously  differentiable  on  its  domain 
of  definition. 

Observe  that  a  thermoelastic  material  is  defined  in  a  manner  such  that  the 
rate  of  entropy  production  i]oc(-;‘P)  in  a  region  V  due  to  mechanical  dissipation 
is  identically  zero  on  T.  Hence,  the  localization  of  the  imbalance  law  (2.1.10) 
at  a  point  contained  in  the  interior  of  M  \  E  yields,  with  the  aid  of  (2.3.2),  the 
inequality 

q(-,  •,  d)  •  d  >  0  on  f+x^  Vd  6  M3  (2.3.3) 

as  a  condition  necessary  for  the  satisfaction  of  the  second  law  of  thermodynamics. 
The  response  function  q  is  assumed  to  be  specified  so  that  (2.3.3)  holds;  then, 
inequality  (2.1.11)$  is  automatically  satisfied  and  can  be  ignored  in  the  following. 

For  remarks  regarding  the  consequences  of  objectivity  on  the  properties  of 
the  potential  tp  and  response  function  q,  see  Jiang  [11]. 

2.4.  Constitutive  specialization.  To  facilitate  the  ensuing  analysis  sup¬ 
pose,  henceforth,  that  the  homogeneous  thermoelastic  body  B  is  thermomechani- 
cally  isotropic.  Then  the  Helmholtz  free  energy  potential  ip  and  heat  flux  response 
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function  q  can  depend  on  the  deformation  gradient  F  only  through  the  deforma¬ 
tion  invariants  J*(G)  defined  in  (2.1.5).  Assume  henceforth  that  both  xp  and  q 
axe  independent  of  the  second  deformation  invariant  /2(G)-  Suppose,  moreover, 
that  the  Helmholtz  free  energy  potential  xp  can  be  represented  in  terms  of  three 
functions  xp  :  2R+  x  1R+  — »  JR,  g  :  JR+  x  2R+  — >  2R  and  g  :  1R+  x  JR+  —*  1R  in  the 
form 

&F,  9)  =  <«!,(  G),  9)  +  *,.(/,  (G), 9)  ft  J3(G), 9)  +  ft/3(G),  9) 

V(F,9)€£+ xffi+,  (2.4.1) 

and  that  the  heat  flux  response  function  q  can  be  expressed  in  terms  of  a  function 
<j> :  JR+  x  JR+  x  JR+  — >  JR  via 

q(F,M)  =  4>(h(G),I3(G),0)d  V(F,0,d)  G  £+  x  x  J23.  (2.4.2) 

In  (2.4.1)  and  (2.4.2)  G  is  the  left  Cauchy-Green  tensor  defined  in  terms  of  the  de¬ 
formation  gradient  tensor  F  by  (2.1.4).  In  accordance  with  the  stipulated  smooth¬ 
ness  of  xp  and  q,  the  functions  xp,  g  and  g  are  taken  to  be  once  continuously 
differentiable  and  piecewise  twice  continuously  differentiable  on  1R+  x  2R+,  while 
<t>  is  taken  to  be  continuous  and  piecewise  twice  continuously  differentiable  on 
1R+  x  1R+  x  1R+.  Assume,  in  addition,  that  the  functions  xp,  g  and  g  comply,  for 
each  6  in  1R+,  with  the  following  isochoric  restrictions: 

*(M)=0,  g(  1,0)  =  0,  $/,(M)  =  - 1,  f(M)  =  J/,(M)  =  0.  (2.4.3) 

In  what  follows,  attention  will  be  restricted  to  homogeneous  isotropic  thermoe¬ 
lastic  materials  wherein  the  Helmholtz  free  energy  potential  xp  obeys  (2.4.3).  A 
particular  material  of  this  type  was  studied  by  Jiang  [11]. 

The  nominal  stress  response  of  B  is  then  determined,  with  the  aid  of  (2.3. 1)2 
and  (2.4.1),  by 


S  =  2^(x.(/i(G),/3(G),9)F  +  Xj(/i(G),A(G),9)F-1-)  on  M\S,  (2.4.4) 
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where  the  functions  \a  JR+  x  1R+  x  2R+  — ► JR  are  defined  as  follows: 

Xi(/i(G),*(G),»)  =  */,(/.,*) +  ^,l/l(/, ,«)«/,,») 

V(/i,/3,d)  6  1R+  x  .R+  x  2R+ , 

(2.4.5) 

V(/i,/3,^)  6  JR+  x  .ffLf  x  JR+. 

Observe  that,  for  an  isochoric  deformation — where  I3  =  1  on  M  \  X1,  use  of  (2.4.3) 
in  (2.4.5)  reduces  (2.4.4)  to  read 

S  =  2pV>j1(Ji,0)  (F  —  F-t)  on  A4\27.  (2.4.6) 

Following  the  work  of  Jiang  &:  Knowles  [12]  in  the  purely  mechanical  set¬ 
ting,  it  can  be  readily  shown  that  a  special  thermoelastic  material  of  the  type 
characterized  by  (2.4.1)-(2.4.3)  satisfies  the  Baker- Ericks en  inequalities  at  all  ab¬ 
solute  temperatures  if  and  only  if 

MIu8)  +  $ilil(hJ)g{h,e)>  0  V(/,j3,e)6WxlR+,  (2.4.7) 

where  the  set  U  is  given  by 

U  =  {(/,,/, )|0  </,<(/, /3)3}. 

Choose  a  rectangular  Cartesian  frame  X  =  {0;ej,e2,e3}  and  consider  the 
response  of  the  thermoelastic  material  at  hand  to  a  simple  shear  deformation  y 
given  as  follows 

y(x, <)  =  (!  +  7e3  ®  ei )x  V(x,t)  €  M,  (2.4.8) 


where  the  constant  7 — assumed  non-negative  without  loss  of  generality — denotes 
the  amount  of  shear.  Note  that  the  foregoing  deformation  is  isochoric.  From 
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(2.3.1) 2,  (2.4.1),  (2.4.3)  and  (2.4.6)  the  nominal  shear  stress  corresponding  to  the 
deformation  y  is,  therefore,  for  each  7  in  M+,  found  to  be 

e3  •  Sei  =  2/ryt/)/, (3  +  72, 9)  =:  r( 7, 9).  (2.4.9) 

where  9  takes  on  some  positive  value.  The  function  r  :  M+  x  JR+  — *  JR  will 
be  referred  to  as  the  shear  stress  response  function  of  the  special  thermoelastic 
material  at  hand  in  simple  shear.  An  immediate  consequence  of  (2.4.3)i  and 

(2.4.9)  is  that  pip  can  be  expressed  via 

v/TT-3 

pj>(h,9)  =  J  T(K,6)dK  V(Ji,0)€  [3,oo)  x.R+,  (2.4.10) 

0 

so  that  the  nominal  stress  response  of  such  a  material,  in  all  three  dimensioned 
deformations  and  absolute  temperatures,  is  completely  characterized  by  specifying 
a  shear  stress  response  function  r  along  with  the  functions  g  and  g  introduced  in 

(2.4.1) .  Now,  define  the  secant  modulus  in  shear  M  :  Bi+  x  JR+  —*  1R  by 

M(7,0)  =  2prjIl(3  +  'r\9)  V(7,0)  G  R+  x  JR+.  (2.4.11) 

Observe  that,  in  compliance  with  the  stipulated  smoothness  of  ip,  both  r  and  M 
must  be  continuous  and  piecewise  continuously  differentiable  on  1R+  x  1R+.  From 

(2.2.9)  and  (2.2.11)  that  the  shear  stress  response  function  r  must  also  satisfy 

r(O,0)  =  0  V0  G  1R+,  r7(0,^)  =  M(0,6)  V«Gfl+.  (2.4.12) 

Note,  also,  that  for  the  simple  shear  deformation  defined  via  (2.1.1)  and  (2.4.5), 
the  Baker- Ericksen  inequality  (2.4.10)  reduces,  with  the  aid  of  (2.4.11)  and 
(2.4.3)4,  to  a  relation  which  involves  only  M:  viz., 


M(7,S)>0  V(7,0)€tf?.. 


(2.4.13) 
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Restrict  attention  in  the  sequel  to  these  special  materials  for  which  the  infinites¬ 
imal  shear  modulus  is  positive;  *.e.,  require  that 

M{ 0,0)  >0  V0  6JR+-  (2.4.14) 

Despite  the  significant  restrictions  which  have  been  placed  upon  the  class 

of  materials  which  will  be  considered  in  this  investigation,  no  presuppositions 

have  been  made  regarding  the  sign  of  the  derivative  with  respect  to  its  first 

argument — where  it  exists — of  the  shear  stress  response  function  corresponding  to 

the  thermoelastic  material  defined  in  compliance  with  (2.4.1)-(2.4.3)  and  (2.4.10). 

Jiang  [11]  has  shown  that  the  monotonicity  of  r(*,  6)  is,  for  fixed  6  in  2R+,  related 

directly  to  the  ellipticity  of  the  material  which  it  characterizes.  If,  in  particular, 

is  not  a  monotonically  increasing  function  on  its  domain  of  definition — for 

some  range  of  9 — then  the  associated  material  is  non-elliptic.  With  this  in  mind, 

let  ( 9  mj^w)  be  contained  in  and  define  functions  “y  •  and 

* 

7  :  (^m,  9m)  M+  such  that 

7(0)  <  7(0)  V0G(0m,0Af).  (2.4.15) 

Next,  define  three  plane  open  subsets  A\,  A2,  and  A3  of  the  shear  strain- 
temperature  quadrant  as  follows: 

=  {(7,0)10  <  7  <  7(0), 0  G  (0m, 0m)}, 

* 

M  =  {(7,0)17(0)  <  7  <  7(0), 0  €  (0m,0M)},  (2.4.16) 

* 

=  {(7, «)!-!(«)  <7  <  °°,0  e 

This  investigation  will  make  use  of  a  particular  subclass  of  non-elliptic  thermoe¬ 
lastic  materialsof  the  above  special  form  wherein  the  relevant  shear  stress  response 
function  r  is  taken  to  be  continuous  on  2R+  x  JR+  and  continuously  differentiable 
on  Ai  U  A2  U  A3  and  is  required  to  obey  the  following  monotonicity  requirements 


r7  >  0  on  Ai  U  A3 , 
r7  <  0  on  A2. 


(2.4.17) 
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Assume,  also,  that  r(*,0)  is  monotonically  increasing  on  for  all  6  in  JR+  \ 
I^mi  &m].  Let  the  nominal  conductivity  in  shear  k  :  Six  ]R+  —*  ]R+  and  the  nom¬ 
inal  specific  heat  per  unit  mass  in  shear  c :  IRx  M+  — ►  1R.+  of  the  thermoelastic 
material  at  hand  be  defined  as  follows 


£(7,0)  =  <£(3  +  72,M)  V(7,0)e  MxR+, 
c(-y,  6)  =  -B4>ee(Z  +  72 , 0)  V(7, 6)  e  JR  x  R+ . 


(2.4.18) 


Suppose  that  k  and  c  are  both  continuous  on  A\  UA2UA3  and  piecewise  continuous 
on  M+  x  IR-f .  In  compliance  with  (2.3.3)  let  k  be  positive  on  its  domain  of 
definition.  Suppose,  in  addition,  that  c  is  positive  on  its  domain  of  definition. 

The  sets  A\  and  A3  are  referred  to  as  the  high  and  low  strain  phases  of  the 
thermoelastic  material  specified  by  (2.4.1)-(2.4.3)  and  (2.4.10).  These,  together 
with  the  set  of  shear  strain-temperature  pairs  in  (M+  x  1 R+)  \  {A\  U  A2  U  A3) 
comprise  the  elliptic  phases  of  such  a  material.  A  thermoelastic  material  of  the 
type  at  hand  which  is  defined  so  that  r,  k  and  c  have  the  properties  set  forth 
above  will  be  referred  to  herein  as  a  three-phase  thermoelastic  material.  See 
Figure  1  for  a  graph  of  t(-,0 0)  for  fixed  90  in  (#m,  #a/)  typical  of  those  which 
specify  three-phase  thermoelastic  materials.  Consult  Figure  2  for  a  division  of 
the  shear  strain-temperature  quadrant  into  regions  of  monotonicity  of  r(-,  6)  for 
fixed  6. 

2.5.  Completion  of  constitutive  assumptions  via  the  kinetic  rela¬ 
tion.  Let  B  be  composed  of  a  three-phase  thermoelastic  material  and  consider 
a  motion  of  B  which  involves  a  moving  surface  of  discontinuity  St  in  one  or  all 
of  the  field  quantities  F(-,<),  u(-,t),  S(-,t),  ^(-,t),  q(-,t),  and  q(*,<)  at  each  in¬ 
stant  t  in  T.  Assume  that  St  separates  high  and  low  strain  elliptic  phases  in  the 
material  at  hand.  In  the  context  of  such  a  motion  it  is  necessary  (see  [1-4])  to 
supplement,  in  some  fashion,  the  constitutive  information  provided  in  Section  2.4. 
An  approach  to  this  taken  by  Abeyaratne  &  Knowles  [1]  entails  the  provision 
of  a  kinetic  relation.  Two  basic  cases  motivated  by  [1]  can  be  considered:  in  the 
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first  a  constitutive  response  function  V  :  JR  x  (0m,0M)  — *■  iR  is  specified  so  that 

Vn  =  V(f,0)  V(f,0)  €  2Rx  (0m,0M),  (2.5.1) 

while,  in  the  second  a  constitutive  response  function  <£  :  JR  x  {6m, 0m)  —*  JR  is 
furnished  so  that 


/  =  0<£(V„,0)  V(V„,0)  eJRx  ( 0m,0M ).  (2.5.2) 

The  functions  F  and  <£  are  referred  to  as  kinetic  response  functions.  Since  the 
three-phase  thermoelastic  material  can  lose  ellipticity  only  for  absolute  tempera¬ 
tures  6  in  (0m, 6m),  the  kinetic  response  functions  V  and  Cp  are  defined  only  on 
JR  x  (0m, 6m)-  Both  varieties  of  kinetic  response  functions  will  be  considered  in 
this  investigation.  If  F  is  such  that  V($,  6)$  >  Ofor  all  (#,  9)  in  JRx(0m,0M  )  then 
(2.2.6)  is  automatically  satisfied  and  F  is  referred  to  as  admissible.  If  <£(F,  6)V  >  0 
for  all  (V,0)  in  JR  x  (0m,0M),  Cp  is,  similarly,  referred  to  as  admissible.  If  an  ad¬ 
missible  kinetic  response  function  V  (or  <£)  is  continuous  on  JR  x  (0m,0M ),  then 
it  must  satisfy  F(O,0)  =  0  (or  <£(0,0)  =  0)  for  all  6  in  (0m,0M)-  If,  in  addition, 
to  being  admissible,  F  (or  Cp)  is  continuously  differentiable  on  JR  x  (6m,  6m),  then 
V#(O,0)  >  0  and  V©(O,0)  =  0  (or  <£y(0, 6)  >  0  and  <£*(0,0)  =  0)  for  all  0  in 
(6m, 6m) — here  V#  and  <£v  refer  to  the  first  partial  derivatives  of  V  and  <£  with 
respect  to  their  first  arguments  while  Ve  and  p$  refer  to  the  first  partial  deriva¬ 
tives  of  V  and  <£  with  respect  to  their  second  arguments.  Otherwise,  admissiblility 
implies  nothing  with  regard  to  the  sign  of  the  derivative  of  a  smooth  kinetic  re¬ 
sponse  function.  All  kinetic  response  functions  considered  herein  are  assumed  to 
be  admissible.  See  Figure  3  and  Figure  4  for  illustrative  graphs  of  V(-,0o)  and 
<£(*,0o)  for  fixed  0o  in  (0m,0Af)- 

2.6.  Thermoelastic  antiplane  shear  motions  of  a  special  thermoe¬ 
lastic  material.  Suppose,  from  now  on,  that  R  is  a  cylindrical  region  and  choose 
a  rectangular  Cartesian  frame  X  =  {0;ei,e2,e3}  so  that  the  unit  base  vector  e3 
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is  parallel  to  the  generatrix  of  71.  A  dynamical  process  will  be  referred  to  as  a 
thermoelastic  antiplane  shear  normal  to  the  plane  spanned  by  the  base  vectors  ei 
and  e2  if  the  deformation  y  is  of  the  form 

y(x,<)  =  x +  «(xi,x2,0e3  V(x,t)  6  M,  (2.6.1) 

and  the  nominal  Helmholtz  free  energy  per  unit  mass  rp,  nominal  entropy  per 
unit  mass  77,  nominal  absolute  temperature  9,  and  nominal  heat  flux  vector  q 
are — like  the  displacement  field  associated  with  (2.6.1) — independent  of  the  x3- 
coordinate.  The  non-trivial  component  of  displacement  u  in  (2.6.1)  will  be  referred 
to  as  the  out-of-plane  displacement  field.  Inspection  of  (2.6.1)  reveals  that  any 
discontinuities  in  the  gradient  and,  perhaps,  time  derivative  of  y  must  result 
from  discontinuities  in  the  spatial  derivatives  out-of-plane  displacement  field  and, 
hence,  occur  across  surfaces  which  do  not  vary  with  the  x3-coordinate;  similarly, 
because  of  their  independence  of  the  x 3-coordinate,  any  discontinuities  in  xj>,  77  or 
q  must  occur  across  surfaces  which  do  not  vary  with  the  x3-coordinate.  Let  St 
denote  a  surface  across  which  at  least  one  of  the  above  field  quantities  jumps  at 
the  instant  t  in  T  and  let  E  be  defined  as  in  (2.1.2). 

Following  the  work  of  Jiang  [11]  in  the  inertia-free  context,  it  is  possible 
to  demonstrate  that,  although  not  every  homogeneous  and  isotropic  thermoelas¬ 
tic  material  can  sustain  thermoelastic  antiplane  shear  motions,  all  thermoelastic 
materials  defined  in  compliance  with  (2.4.1)-(2.4.3)  and  (2.4.10)  are  capable  of 
doing  so.  It  is  easily  shown  that  for  such  materials  the  local  balance  equations 
(2.1.11)2,3,4,5  reduce,  in  the  absence  of  body  forces  and  heat  supplies,  to 

(M(7, 6)u,a ) ,«  =  pn  on  X\T, 

(2.6.2) 

{k(7,B)9,a),Q+M9{~t,0)6u,aU,a=  pc(l,9)9  on  X\T, 

where  X  is  given  by  V  x  T,  V  is  a  generic  cross  section  of  71,  and  F  = 
{(x1,x2,<)|(xi,x2)  €  Ct,t  €  T]  with  Ct  =  T>f)St  at  each  t  in  T.  See  Fosdick  & 
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Serrin  [8]  and  Fosdick  &  Kao  [7]  for  a  general  discussion  of  the  circumstances 
under  which  the  local  balance  equations  (2.1.11)2,3  reduce  to  a  single  scalar  equa¬ 
tion.  In  (2.4.2)  M  is  the  secant  modulus  in  shear  as  defined  in  (2.4.11)  and 
7  :  X  \  r  — >  M  is  the  shear  strain  field  given  by 

y(xux2,t)  =  \/u,a(xi,X2,t)u,a(xi,X2,t)  V(xi  ,x2,t)€X\r.  (2.6.3) 

The  jump  conditions  (2.1.13)2,3  reduce,  for  a  thermoelastic  material  of  the 
type  at  hand  subjected  to  antiplane  shear,  to 


W( 7, 0)u,a  na]  +  pvn[u\  =  0  on  J1, 
,0)0, a  »*o]  +  -  fVn  =0  on  r, 


(2.6.4) 


where  F  —  {(x,t)|x  6  Ct,t  €  T),  n(*,t)  :  Ct  — *  M  is  a  unit  normal  to  Ct,  the 
nominal  entropy  per  unit  mass  T)  :  X  \  T  —*  Mis  given,  from  (2.3. 1)3,  (2.4.1)  and 
(2.4.3),  by 


7 

n  =  -^«(3  +  72,tf)  =  J  Te(K,0)dK  on  X  \ 


(2.6.5) 


and  f  :  r  — ►  1R  is  the  driving  traction  introduced  in  Section  2.3.  The  kinematic 
jump  condition  (2.1.16)  becomes 


fu]  =  0  on  r,  [^J  =  0  on  J1. 


(2.6.6) 


It  is  also  readily  shown  that  the  driving  traction  /  for  a  thermoelastic  material 
defined  via  (2.4.1)-(2.4.3)  and  (2.4.10)  subjected  to  an  antiplane  shear  deforma¬ 
tion  involving  a  discontinuity  in  the  gradient  of  the  out-of-plane  displacement 
across  a  moving  curve  Ct  is  given  by 


+ 

=  J_  r(K,9)dK  -  ((M(y,  0)u,a))[u,o] 


on  r. 


(2.6.7) 
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Recall  that  the  jump  condition  (2.1.13)4,  or,  equivalently,  (2.2.7)  is  satisfied  con- 
stitutively  by  requiring  that  the  kinetic  response  function  be  admissible. 

With  reference  to  (2.1.14),  (2.1.15)  and  (2.2.6)  it  is  easily  demonstrated  that, 
in  the  absence  of  inertial  effects,  (2.4.2)  is  replaced  by 

(M(7,0)u,q),q  =  O  on  X\T,  (2.6.8) 


while  (2.4.4)  becomes 


\M( 7, 8)u,a  na]  =  0  on  T, 


(2.6.9) 


and  (2.6.7)  reduces  to 

+ 

/  =  /  t(k,6)cIk  -  M(7,6)u,q  [u,a]  on  T.  (2.6.10) 

Ji 

Observe  that,  within  the  context  of  a  thermoelastic  antiplane  shear  defor¬ 
mation  of  the  type  described  above,  no  generality  is  lost  by  focussing  exclusively 
upon  the  motion  on  a  cross-section  T>  of  the  cylinder  7 Z  and  the  dynamics  of  the 
curve  Ct  =  V  fl  St-  In  the  following,  curves  Ct  across  which  the  gradient  of  the 
out-of-plane  displacement  field  u(-,-,t)  and,  perhaps,  the  out-of-plane  velocity 
field  «(-,-,<),  the  entropy  field  !](•,-,  t),  and  the  gradient  of  the  absolute  temper¬ 
ature  field  0(‘,',t)  jumps,  at  some  instant  t  in  T,  and  which  segregate  the  high 
and  low  strain  phases  of  the  material  at  hand  will,  therefore,  be  referred  to  as 
phase  boundaries. 
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3.  LINEAR  STABILITY  OF  A  PROCESS  INVOLVING  A  STEADILY 
MOVING  PLANAR  PHASE  BOUNDARY  IN  A  THREE-PHASE 
THERMOELASTIC  MATERIAL 

3.1.  Description  of  the  base  process.  Suppose  that  B  is  composed  of 
a  three-phase  thermoelastic  material  and  that  the  cylinder  %  degenerates  so  as 
to  occupy  all  of  2R3.  Let  the  rectangular  Cartesian  frame  X  be  as  in  Section  2.4. 
Consider  a  thermoelastic  antiplane  shear  motion  on  the  time  interval  (— oo,0) 
with  an  out-of-plane  displacement  field  uo(-,t)  :  1R  — >  1R  given  by 


7jii  +  V[t 
JrXl  +  M 


if 

if 


Xl  <  Vot, 
Xi  >  Vo t, 


(3.1.1) 


for  each  t  in  (— oo,0),  and  an  absolute  temperature  field  6q  which  is  constant  on 
It  x  (—oo,0)  and  satisfies 

Ooe(0m,OM),  (3.1.2) 

where  the  shear  strain-temperature  pairs  (71,60)  and  (7,-,  60)  satisfy  one  of  the 
following 


((7/1^0), (7r,$o))  €  A3  x  Ai,  ((7ij0o),(jr,0o))  €  A 1  x  A3.  (3.1.3) 

Observe  that  the  process  described  by  (3.1.1)— (3.1.3)  is  isothermal. 

Since  one  of  (3.1.3)  must  hold,  there  is  no  loss  in  generality  incurred  by 
assuming  that  the  base  interface  normal  velocity  vq  is  non-negative;  that  is, 


v0  >  0.  (3.1.4) 

It  is  clear  that  uq  and  9q  satisfy  the  differential  equations  in  (2.6.2)  on  the 
set  (2R2  x  (-oo,0))  \  T0  with  JT0  given  by  {(xi,x2,0l  (^1^2)  €  At,t  6  (-00, 0)} 
and  At  =  {(2:1,12)!  ^1  =  vot,X2  €  1R.}  for  each  t  in  (— 00, 0).  The  moving  line  At 
is,  for  each  t  in  (—00,  o),  a  phase  boundary. 
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Assume,  in  order  to  comply  with  the  jump  conditions  in  (2.6.4)  and  (2.6.6) 
on  To,  that  the  constants  7 /,  7r,  v/,  vr,  and  vo  associated  with  (3.1.1)— (3.1.3)  are 
restricted  to  satisfy  the  following  equations: 

Vr  -  VI  -f  U0(7r  -  7 0  =  0, 

T(7r,0o)  -  t(7j,  0o)  +  pvo(vr  ~  Vj)  =  0,  (3.1.5) 

Voifo  +  p0o(t]r  -  Vl))  =  0. 

In  (3.1.5)  the  base  driving  traction  /o  is  given,  with  the  aid  of  (2.6.7),  by 

“Yr 

fo  =  /r(7,,o)d7-i(.(7^o)  +  r(7^o))(7r-70,  (3.1.6) 

■u 

and  the  constants  ?/r  and  rji  are  given  in  terms  of  7r,  7/  and  do  via  the  shear  stress 
response  function  r  as  follows 

It  'll 

Vr  =  --  J  Tg(K,d0)dK,  iy  =  J  Te(K,d0)d.K.  (3.1.7) 

0  0 

Observe,  as  a  consequence  of  (3.1.2)  and  (2.2.7),  that  fo  must  satisfy 

/o  >  0.  (3.1.8) 


Assume  that  vq  complies  with  the  inequality 


vo  <  min 


(3.1.9) 


so  that  the  normal  velocity  of  the  phase  boundary  in  the  base  process  is  locally 
subsonic.  It  is  then  permissible2  to  impose  a  kinetic  relation  in  the  form  (2.3.8) 


2  See  Abeyaratne  &  Kknowles  [2]. 
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or  (2.3.9)  on  Jo  and  require  that  the  parameters  7 /,  7 r,  vj,  vr,  and  Vo  satisfy  one 
of 


vo  =  V(£,0o),  /o  =  Qo<p(vo,0o),  (3.1.10) 


depending,  respectively,  upon  whether  a  kinetic  relation  is  provided  in  the  form 
(2.5.1)  or  (2.5.2). 

In  a  coordinate  frame  moving  with  the  phase  boundary,  the  base  process 
described  involves  a  piecewise  homogeneous  shear  strain  field  and  a  homogeneous 
temperature  field.  If  (7/,  So)  and  (7,-,  do)  are  consistent  with  (3.1.3)i  then  the  base 
process  is  one  wherein  the  high  strain  elliptic  phase  of  the  material  at  hand  grows 
at  the  expense  of  the  low  strain  elliptic  phase  at  constant  temperature;  whereas,  if 
(7/,  60)  and  (7r, 60)  comply  with  (3. 1.3)2  then  the  base  process  is  such  that  the  low 
strain  elliptic  phase  of  the  material  at  hand  grows  at  the  expense  of  the  high  strain 
elliptic  phase  at  constant  temperature.  In  either  case  the  discontinuity  involved 
is,  for  the  duration  of  the  motion,  a  normal  phase  boundary — that  is,  the  angle 
between  the  limiting  values  of  the  gradient  of  the  out-of-plane  displacement  field 
on  either  side  of  the  phase  boundary  is  zero  at  every  point  of  the  phase  boundary 
over  the  time  interval  (-oo,0). 

The  constant  latent  heat  of  transformation — Co — associated  with  the  ther¬ 
moelastic  process  described  by  (3.1.1)— (3.1.3)  is  defined  by 


Co  =  pdoim  -  V r)  -  fo  =  Oo 


'll 


(3.1.11) 


From  (3.1.5)3  it  is  clear  that  C0  must  be  zero  if  vo  =  0 — which  agrees  with  the 
intuitive  notion  that  the  heat  given  off  in  the  process  of  transformation  must 
be  zero  in  the  absence  of  heat  flux.  Recall  from  Section  2.5  that,  under  the 
present  assumption  that  the  kinetic  response  function  V  or  ip  which  is  provided 
is  continuous,  Vo  =  0  if  and  only  if  fo  =  0.  Hence,  when  vq  =  0,  the  latent  heat 
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of  transformation  simplifies  to 

7r 

to  =  p9o(rn  —  r]r)  =  80  J  Te(K,6o)dK.  (3.1.12) 

71 

Observe,  however,  that  (3. 1.5)3  is  satisfied  for  any  real  value  of  to  when  vo  =  0. 

Suppose,  in  addition  to  all  the  above,  that  the  kinetic  response  function  V 
or  <p  is  chosen  so  that  its  derivative  is  non-zero  at  the  base  driving  traction  /o, 
that  is,  assume  that  one  of  the  following,  as  is  appropriate  to  the  specification  of 
a  kinetic  relation  in  the  form  of  either  (2.5.1)  or  (2.5.2),  must  hold: 

v*(£ A)^o,  <Mv0A)/o,  (3.1.13) 

This  assumption  is  made  in  order  to  preclude  the  necessity  of  going  to  higher 
order  in  the  context  of  the  forthcoming  linear  stability  analysis.  See  Figure  3  and 
Figure  4  for  schematic  graphs  of  smooth  admissible  kinetic  response  functions 
V(-,0q)  and  <p(-,0o)  which  satisfy  (3.1.13). 

When  inertial  effects  are  ignored  it  is  clear  that  uo  as  defined  in  (3.1.1)  also 
satisfies  the  field  equation  in  (2.6.9)  on  ( 1R 2  X  (— oo,0))  \i"b.  Equations  (3.1.5)i,3 
are,  in  this  context,  still  sufficient  to  satisfy  (2.6.8)  and  (2.6. 2)2  on  JT0-  In  place 
of  (3. 1.5)2,  the  constants  7 /,  7r,  V(,  vr,  and  Vo  must,  however,  satisfy 

T(lr,  So)  -  r(7f,  0O)  =  0,  (3.1.14) 

in  order  for  the  jump  condition  in  (2.6.9)  to  hold  on  Jo-  Although  the  expression 
for  the  base  driving  traction  fo  given  in  (3.1.6)  remains  valid  in  the  inertia-free 
setting,  (3.1.14)  can  in  this  case  be  used  so  that  it  simplifies  to  read 

fo  =  J  r( 7,  S0)  dj  -  r.( 7r  -  7 1), 

71 


(3.1.15) 


where  r.  =  r(7/,0o)  =  T(yr,0o). 

Given  a  shear  stress  response  function  r  which  describes  a  particular  three- 
phase  thermoelastic  material  and  an  arbitrary  kinetic  response  function  V  or 
Cp  which  describes  the  dynamics  of  phase  boundaries  which  may  occur  therein, 
there  may  not,  in  general,  exist  constants  7/,  7r,  vi,  t>r,  and  i>o  which  satisfy  one  of 
(3.1.3)i  or  (3. 1.3)2  and  are  consistent  with  the  restrictions  embodied  by  (3.1.5), 
(3.1.9),  (3.1.10)i  or  (3.1.10)2,  and  (3.1.13)i  or  (3.1.13)2,  or,  in  the  inertia-free 
case,  (3.1.5)i,  (3.1.13)i  or  (3.1.13)2,  (3.1.9),  (3.1.10)i  or  (3.1.10)2,  (3.1.14)  and 
(3.1.15).  Within  the  context  of  this  investigation  it  will  be  assumed,  however, 
that  V  or  <p  is  chosen  so  that  a  non-trivial  base  process  exists. 

3.2.  Perturbation  of  the  base  process.  Suppose  that  at  the  instant  t  =  0 
the  out-of-plane  displacement  and  velocity  fields,  the  absolute  temperature  field 
and  the  configuration  of  the  phase  boundary  associated  with  the  thermoelastic 
process  specified  in  Section  3.1  are  subjected  to  a  perturbation.  Let  this  perturba¬ 
tion  be  such  that  the  phase  boundary  can  be,  at  t  =  0+,  described  by  the  graph 
Co  of  a  continuous  function  h  :  JR  — ►  JR  of  the  a^-coordinate,  and  segregates 
elliptic  phases  of  the  three-phase  material  at  hand  in  a  sense  consistent  with  that 
which  was  present  for  t  in  (— 00, 0).  Let  the  out-of-plane  displacement  field,  out- 
of-plane  velocity  field,  and  absolute  temperature  field  linked  to  this  perturbation 
be  given,  respectively,  by  a  once  continuously  differentiable  function  t)  :  JR2  —*  JR, 
a  continuous  function  va  :  JR.2  — >  ]R,  and  a  continuous  function  <f>  :  JR 2  —>  JR. 
Assume  that  h,  rj,  w  and  <j>  represent  small  deviations,  in  some  appropriate  sense, 
from  their  counterparts  in  the  base  process.  In  particular,  suppose  that  h,  rj, 
r),Q,  w,  and  4>  are  all  square  integrable  on  their  domains  of  definition.  Require, 
furthermore,  that  the  components  of  the  gradient  of  7/  allow  the  satisfaction  of 
the  decay  condition 
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while  a?  and  4>  comply  with  the  following  decay  conditions 

lim  w(xi,x2)  =  0,  lim  <£(xi,x2)  =  0,  (3.2.2) 

oo  **+*?— oo 

so  that  the  disturbance  is  localized  in  a  neighborhood  of  the  phase  boundary 
associated  with  the  base  state  at  t  =  0. 

The  perturbation  at  t  =  0  will  initiate  a  new  process  involving  an  out-of- 
plane  displacement  field  u  :  JB?  x  JR+  — »  1R  and  an  absolute  temperature  field 
9  :  IF?  x  M+  — ►  IR  which  are,  in  general,  weak  solutions  of  the  field  equations 
(2.6.2)  and  satisfy  the  jump  conditions  in  (2.6.4)  and  (2.6.6)  at  all  discontinuities 
in  their  gradients,  the  kinetic  relation  (2.5.1)  or  (2.5.2)  at  all  phase  boundaries, 
and  the  initial  conditions 


it(-,-,0+)  =  u0(-,0+)  +  T]  on  M2, 

u(‘,  •,  0+)  =  tto(‘,  0+)  +  w  on  JR2 ,  (3.2.3) 

fl(-,-,0+)  =  9o  +  $  on  R2. 

Since  the  perturbation  is  small,  assume  that,  the  subsequent  process  involves 
only  a  single  phase  boundary  Ct  =  {(xi,X2,t)|xi  =c(x2,<), X2  6  R}  for  each  t  in 
R+ ,  with  q  :  Rx  R+  —*  R  continuously  differentiable  on  its  domain  of  definition 
and  defined  so  that  it  is  in  accord  with  the  initial  condition 

?(-,0+)  =  h  on  R.  (3.2.4) 

With  the  intent  of  linearizing  the  field  equations  in  (2.6.2)  about  the  base  process, 
write,  for  each  t  in  JR+, 

u(xj ,  x2,  t)  =  u0(xi,t)  +  w(*i,x2  ,<)  V(xi,x2)  eT>\Ct , 

(3.2.5) 

9(x i ,  x2,  t)  =  So  +  T(xj ,  X2,  t )  V(xi ,  X2 )  €V\Ct, 
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where  w  and  its  derivatives  and  T  are  assumed  to  represent  small  departures  from 
the  relevent  quantities  in  the  base  process.  Assume  that  the  components  of  the 
gradient  of  w  satisfy  the  following  limits 

lim  tu,i  (xi,-,*)  =  0  on  1R  x  #l+, 

ii— *±oo 

lim  w,2  (',X2>')  =  0  on  JRxM+, 

x  3-*±oo 

and  also  that  T  complies  with  the  limit 

lim  T(xi,X2,-)  =  0  on  2Ft+. 

**+**— »oo 

From  (3.2.3)  and  (3.2.5)  it  is  clear,  moreover,  that — when  inertial  effects  are  not 
ignored — the  increment  w  to  the  out-of-plane  displacement  field  must  satisfy  the 
following  initial  conditions: 


(3.2.6) 


(3.2.7) 


u>(-,  -,0+)  =  V  on  1R2, 
u>(v,0+)  =  G7  on  ]R2. 


(3.2.8) 


It  is  important  to  emphasize  that  these  can  not  be  imposed  in  the  inertia-free 
setting. 

Also,  the  increment  T  to  the  absolute  temperature  field  must  satisfy  the 
following  initial  condition 


T(-,-,0+)  =  <t>  on  JR2. 


(3.2.9) 


Next,  define  s  :  JF2  x  .51+  — ►  1R,  the  correction  to  the  interface  position  due 
to  the  perturbation,  via 


=  vot  +  on  Si  V<  €  -5£+. 


(3.2.10) 
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Note,  from  (3.2.4)  that  the  increment  s  to  the  phase  boundary  position  must 
satisfy  the  initial  condition 


s(-,0+)  s=  h  on  1R.  (3.2.11) 

Observe  that  the  unit  normal  vectors  n±(*,  t) :  2R  —*  M  to  Ct  are  given  by 
D±(-,f)  =  ±6’ 3,2  asaassr  Oil  M  Vf  €  R+.  (3.2.12) 

Vi +  «,!(•,() 

For  the  remainder  of  this  work,  choose  the  unit  normal  vector  associated  with  the 
plus  sign  in  (3.2.12)  and  drop  this  sign  when  referring  to  it.  The  normal  velocity 
Vn(-,<)  :  ]R  —¥  M  of  Ct  is,  thus,  given,  for  each  t  in  .R+,  by 

V„(-,t)  =  ?  +  K  V(€B+.  (3.2.13) 

V 1  +  si2  (‘>0 

3.3.  Linearization  of  the  field  equations  associated  with  the  process 
initiated  by  the  perturbation.  Let  T>[  and  V\  denote,  for  each  t  in  2R+,  plane 
sets  defined  as  shown  below: 

Vlt  =  {(xi,x2)\x1  <<(x2,f)},  v;^m2\i>i.  (3.3.1) 

Let  Xi  and  Xr  be  given,  in  turn,  by 

Xi  =  {(xi,X2,<)|(xi,X2)  €  V[,t  €  1R+),  (3.3.2) 

and 

Xr  =  {(xi,X2,<)|(xi,X2)  G  G  1R+}.  (3.3.3) 

O  O 

The  displacement  equations  of  motion  which  hold  on  Xi  and  Xr  can  be  obtained 
following  Fried  [9]  and  are  given,  in  turn,  by 

a?tu,  u  +tfw,  22  = 

a2rw,u  +b\w,22  =  ti>. 


(3.3.4) 
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where  the  positive  constants  a/,  6j,  and  ar,  br  are  defined  as  follows: 

at  =  \/Ty('ri,Oo)/p,  bi  =  y/M(ji,60)/p , 

V _  _  (3-3.5) 

<*r  =  y^y(7rA)/p,  =  y/M(jr,0O)/ p. 

In  writing  (3.3.5),  the  positivity  of  t7(7j,0o)  and  r7( 7r,0o) — which  are  results 
of  whichever  of  (3.1.3)i,2  is  appropriate,  and  of  M(yi,  do)  and  M(-yr,0o) — which 
follow  from  (2.4.13),  have  been  used. 

O  O 

The  energy  equations  which  hold  on  Xi  and  Xr  can  be  obtained  by  linearizing 
the  partial  differential  equation  (2. 6. 2)2  about  (7/,  0o)  and  (7r,  0o),  respectively. 

O 

Turn,  now,  to  the  derivation  of  the  linearized  energy  equation  which  holds  on  Xi . 
It  is  easy  to  show,  following  [9],  that 

7  =  7/  +  u>,i  on  Xi.  (3.3.6). 

With  the  aid  of  (3.2.5)2  and  Taylor’s  theorem  the  relation  (3.3.6)  leads  to  the 
following  expansions: 

fc(7,0)  £  fc(7/,0o)  +  k1(yi,60)w,i  +kg(lfl,0o)T  on  X(, 

c('Y,6)  =  c('yi16o)  +  c1(^i,0o)iv,1+cg(ji,0o)T  on  Xt,  (3.3.7) 

Mg( 7, 0)  9£  Mg(yi,  0o )  +  0o)u>,i  +Mw( 7 j,  0o)T  on  Xi. 

Next,  using  (3.2.5)  and  (3.3.7)i,2  in  the  left  hand  side  of  the  partial  differential 
equation  in  (2.6.2)2  gives 

Ck('Y,O)0,a),a  +  M$(y10)0u,ou,a* 

k(7i,0o)T,Qa+0oT9(n,0o)w,i  on  Xt,  (3.3.8) 

while  using  (3.2.5)  and  (3.2.7)3  in  the  right-hand-side  of  the  same  equation  gives 

pc(7, 0)0  3  pc{‘ji,0o)T  on  X\.  (3.3.9) 
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O 

The  linearized  energy  equation  which  holds  on  Xi  is,  thus,  from  (3.2.5),  (3.3.8) 
and  (3.3.9),  given  by 

<xiT,aa  =  T  +  ,  (3.3.10) 

where  the  positive  constant  a/  and  the  real  constant  /3/  are  defined  by 


a  _  0oTg(7t>flo) 
pc(7/,0o)’  1  pc(yt,6Q) 


(3.3.11) 


Similarly,  the  linearized  energy  equation  which  holds 


T  +  0rw, i  » 


(3.3.12) 


where  the  positive  constant  ar  and  the  real  constant  fiT  are  defined  by 


fc(7r,#o)  g  8oTe(lr,0o) 

pc(yr,6oY  r  pc(yr,9o)  ' 


(3.3.13) 


Prom  (2.6.9)  it  is  clear  that,  in  the  inertia-free  setting,  the  displacement 
equations  of  motion  (3.3.4)  are  supplanted  by 


ajjw,u  +tfw,  22  =  0, 
a2rw,\\  +b^w,  22  =  0, 


(3.3.14) 


o  o 

which  hold,  repectively,  on  X\  and  Xr. 

3.4.  Linearization  of  the  jump  conditions  and  kinetic  relation  as¬ 
sociated  with  the  process  initiated  by  the  perturbation.  Since  the  set 
r  =  {(zj,Z2iOI  (xi>x2)  €  Ct,t  €  1R+}  represents  the  post-disturbance  trajec¬ 
tory  of  the  phase  boundary,  the  jump  conditions  in  (2.6.4)  and  (2.6.6)  and  the 
kinetic  relation  in  (2.5.1)  or  (2.5.2) — with  Vn  and  /  given,  respectively,  by  (3.2.13) 
and  (2.6.7) — must  hold  on  it.  Assume,  henceforth,  that  the  function  s  introduced 
via  (3.2.10)  and  its  derivatives  are  small  in  the  same  sense  that  w  and  T  are  small. 
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Note,  first,  that  this  assumption  implies,  using  (3.2.12)  and  (3.2.13),  the  following 
approximations  for  n  and  Vn  on  J1 : 


n  =  ei  —  s,  2  e2  on  F, 
Vn  =  uo  +  s  on  r. 


(3.4.1) 


It  is  easy  to  show,  following  Fried  [9],  that  the  linearized  form  of  the  jump 
condition  (2.6.4)i  is  as  follows 


(' Or  -  vl)w,i  (v0t+,X2,t)  -  (of  -  Vq)w,i  ( V0t~,X2,t ) 

=  2u0(7r  -  n)s(x2,t)  V(x2,<)  €  M  x  2R+, 


(3-4.2) 


while  that  of  (2.6.6)i  is 


w(v0t+,x2,t)  —  w(v0t—,x2,t)  =  (7/  —  'Yr)s(x2,t)  V(x2,t)  €  IRx  M+.  (3.4.3) 

Linearization  of  the  jump  condition  (2. 6. 6)2  gives,  next,  since  no  heat  flux  is 
present  in  the  base  process  (3.1.1)-(3.1.3), 


T(vot-h,  x2,  t)  —  T(uo<— ,  x2,  t)  =  0  V(x2,  t)  6  -R  x  JR+,  (3.4.4) 

so  that  the  increment  T  to  the  absolute  temperature  field  is  continuous  across  the 
phase  boundary  in  the  post  perturbation  process. 

The  driving  traction  /  can  be  linearized  in  a  manner  analagous  to  that  dis¬ 
played  in  [9]  to  give,  with  the  aid  of  (3.4.4), 

/(x2,t)  ^  f0  +  71  -  7r)((a^  -  vl)lV,X  (v0t+,X2,t)  +  (af  -  Uq)^,1  (t>0 <-,X2,t)) 

+  (p(m  -  Vr)  +  5(7 /  -  7r)(T*(7r,  00  )  +  Te^l,  0O)))  T(v0t,  X2,  t ) 

V(x2,t)  €  2R  x  1R+,  (3.4.5) 

where  the  base  driving  traction  /o  is  given  by  (3.1.6). 
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Tum,  now,  to  the  linearization  of  the  energy  jump  condition  (2. 6.4)2.  From 
(3. 3. 7)i,  (3.4.1)i  and  (3.2. 5)2  it  is  clear  that  the  first  term  on  the  left-hand-side 
of  the  energy  balance  jump  condition  (2. 6.4)2  linearizes  as  follows: 

[£(7(C(*2,  t),  X2,  t),d(t(x2,t),  X2,t))8,a  (C(*2 ,  *),  *2, 0n*(*2, 01 

-  fc(7r,0o)7\i  (u0<+,*2,<)  -  £(7/A)T,i  (vQt-,X2,t) 

V(x2,<)  G  1R  X  JR+.  (3.4.6) 

Furthermore,  (2.6.5),  (2. 6.6)2,  (3.4.1),  (3.4.4)  and  (3.4.5)  yield  the  following  lin¬ 
earization  of  the  remaining  two  terms  on  the  left-hand-side  of  the  energy  balance 
jump  condition  (2.6.4)2: 

(p0(((x2,t),  X2,  t)[r/(((X2,t),  X2,t)]  +  f(x2,t))V„(x2,t) 

=pv0(c(7r,  0O)  -  c(y/,0o))T(vot,X2,t) 

+  jMli  -  7r)(Te(yr,0o)  +  re(yl,0o)))T(vot,X2,t) 

+  ipM'TI  -  7rX(«r  -  V%)w,l  (v0t+,X2,t)  +  («?  -  Vo)«M  (v0t-,X2,t)) 

-  vo0o{Tg('yr,0o)w,i  (v0t+,X2,t)  -  T®(7/,0O)^,1  (u0t-,a:2,t)) 

-t0s(x2,t)  V(x2  ,t)£]RxlR+.  (3.4.7) 

Together,  (3.4.6)  and  (3.4.7)  give  the  following  expression  for  the  linearization  of 
the  energy  jump  condition  (2. 6.4)2: 

0  =fc(7r,0o)T,i  (u0t+,Z2,<)  -  fc(7j,0o)7\i  (v0t-,X2,t) 

+  pvo(c(jr,0o)  -  c(jl,0o))T(vOt,X2,t) 

+  $V0(7I  -  7r)(M7r,0o)  +  M7l,0o)))T(vot,X2,t) 

+  \pva(7l  -  7r)((Or  -  Vo)u>,l  {v0t+,X2,t)  +  (dj  -  V%)w,i  (v0t-,X2,t)) 

~  VO0o(To(7r,0o)w,l  (vot+,X2,t)  -  T*(7/,  80)w,i  (l?0 <-,X2,<)) 

—  £o^(*c2,t)  V(x2 ,t)  G  JR  x  2R+.  (3.4.8) 
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Linearization  of  the  kinetic  relation  and  use  of  whichever  of  (3.1.10)  is  ap¬ 
propriate  in  a  manner  completely  analogous  to  that  performed  in  [9]  gives,  with 
the  aid  of  (3.4.4), 


s(x2,t)  — -((a2r  -  vl)w,i  (vQt+,x2,t)  +  (a]  -  i>o)iu,i  ( v0t-,x2,t )) 

+  \  (M^r,e0)  +  Te(n,0o)))T(vot,x2,t) 

lpvm 

+  (v#  +  Pv?e0  )TM,  x2,t)  V(x2,t)  €  2R  x  JR+,  (3.4.9) 


where  the  constants  vm  and  v#  are  defined  by  either 


(3.4.10) 


if  the  kinetic  relation  is  furnished  in  the  form  (2.5.1),  or 


0o<pv(vo,#o)  <pe(vo,Oo) 

p  V>v{vo,8o) 


(3.4.11) 


if  the  kinetic  relation  is  supplied  in  the  form  (2.5.2).  Note,  from  (3.1.13),  that  u* 
is  a  real — but  nonzero — constant,  while  v#  is  a  real — and  possibly  zero — constant. 

By  virtue  of  the  foregoing  calculations  it  is  crucial  to  note  that,  within  the 
scope  of  the  linearization,  it  is  legitimate  to  enforce  the  partial  differential  equa¬ 
tions  in  (3.3.4)i  and  (3.3.10)  on  the  interiors  of  the  set  f2i  defined  by 


/?/  =  {(xi,x2>t)l(*i>*2)  €  n\,t  €  .R+},  (3.4.12) 

with  llj  =  {(xi,X2)ixi  <  vat,x 2  G  JR}  for  each  t  in  1R+,  instead  of  the  set  Xi, 
and  the  partial  differential  equations  in  (3. 3. 4)2  and  (3.3.12)  on  the  interior  of  the 
set  Qr  defined  by 


Qr  =  {(xi,X2,t)|(xi,x2)  e  ni,t  €  -R+}, 


(3.4.13) 
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with  lit  —  {(^i>®2)|xi  >  v0t,x 2  €  2R)  for  each  t  in  1R+,  instead  of  the  set  XT- 
For  the  purposes  of  the  forthcoming  analysis  it  is  useful  to  define  a  set  /  as  follows: 

/  =  {(xi,x2,t)\xi  =  v0t,x 2  €  ]R,t  €  2R+}.  (3.4.14) 

In  the  inertia-free  case  it  is  readily  shown  that,  while  (3.4.3)  and  (3.4.4) 
continue  to  hold,  (3.4.2)  is  replaced  by 

(vot+,x2,t)  —  a;iu,i  (vot—,X2,t)  =  0  V(x2,f)  €  12  x  2R+,  (3.4.15) 

(3.4.8)  reduces  to 

0  =fc(7r,0o)T,i  (v0t+,x2,t)  -  £(7t,0o)T,i  (u0t-,x2,<) 

+  /9Uo(c(7r,^o)  -  c('yi,0o))T(vOt,X2,t) 

+  \v0{~1l  -7r)(r<»(7r,^o)+  Mjl,90)))T(v0t,X2,t) 

+  \pvo{h  -  lr){a2rw,-i  {v0t+,x2,t)  +  a?u;,i  (v0<-,x2,<)) 
-uo6o(re(7r,0o)^,i  (u0*+,x2,<)  -  T*(7i,0o)uM  (vot-,x2,t)) 
-4i(x2,<)  V(x2,t)  €  x  2R+,  (3.4.16) 

and  (3.4.9)  simplifies  to  read 

s(x2,f)  =  {a2rw,i  (v0t+,x2,t)  +  a2w,\  (v0t-,x2,t)) 

+  -ir—^L(Te('fr,So)  +  Tg('yi,9o)))T(v0t,X2,t) 
zpv* 

+  (va  +  "j^0o )T(v°t,x2it)  V(x2,f)  G  JR  x  2R+.  (3.4.17) 

Finally,  remarks  analogous  to  those  made  regarding  the  enforcement  of  the  partial 
differential  equations  in  (3.3.4)i  and  (3.3.4)2  on  hi  and  hr  apply  also  to  those  in 
(3.3.14)x  and  (3.3.14)2. 
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3.5.  Specialization  of  the  base  process  and  the  associated  linearized 
description  of  the  post  perturbation  process.  Suppose,  henceforth,  that  the 
base  process  described  in  Section  3.1  is  restricted  so  that 

re(yi,0o)  =  Mir,  60)  =  0;  (3.5.1) 

it  is  implicitly  assumed  that  the  shear  stress  response  function  r  allows  for  the 
possibility  of  (3.5.1).  Observe  that  (3.5.1)  requires  that  the  coefficients  of  ther¬ 
moelastic  coupling  in  the  low  and  high  strain  phases  of  the  material  at  hand  axe, 
by  (3.3.11)2  and  (3.3.13)2,  both  identically  zero.  Although  this  assumption  is 
made  in  order  to  simplify  the  forthcoming  analysis,  it  is  not  inconsistent  with 
the  isochoric  nature  of  the  deformation  under  consideration.  The  linearized  field 
equations,  jump  conditions,  kinetic  relation,  initial  conditions  (where  appropri¬ 
ate),  and  fax  field  decay  conditions  satisfied  by  the  increments  w,  T  and  s  to 
the  out-of-plane  displacement  field,  absolute  temperature  field  and  the  interface 
position  axe  now  listed  in  both  the  inertial  and  inertia-free  cases. 

In  the  inertial  case,  (3.3.4),  (3.3.10)  and  (3.3.12)  give  the  following  linearized 
field  equations 


a]w,u  +tfw,22  =  w 

on 

tii. 

a2rw,i  1  +blw,  22  =  w 

on 

tin 

otiT,aa  =  T 

on 

tii , 

arT,a(3i  —  T 

on 

tin 

In  addition,  from  (3.4.2),  (3.4.8),  (3.4.3)  and  (3.4.4)  the  following  jump  conditions 
hold 

[(a2  -  Vq)um  1  =  2v0(7r  -  7/)«  on  /, 

V*T,  1  J  +  pv0{clT  =  pvo(ir  -  7/){{(a2  -  Uo)«M  ))  +  *0  *  on  /, 

(3.5.3) 

M  =  (7/  -  7»>  °n  I, 

[0]  =  0  on  /, 
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where  the  constants  a+  and  a2_  are  given  by 


2  2  2  2 
°+  =  ari 


and  k+  =  kr,  k_  =  fc/,  c+  =  cr,  and  c_  =  c/  are  defined  via 


(3.5.4) 


k+  =  fc(7r,0o),  k_  =  k(yi,60), 

C+  =  C(  C_  =  c(7/,^0). 


(3.5.5) 


Next,  from  (3.4.9)  and  (3.5.1)  the  following  linearized  kinetic  relation  holds: 

A  'll  ~  Hz  /// \\  ,  .  l0  \t r  to  K 


^  =  J^<<(«2-^oK1))  +  (^  +  ^)T  on  /. 


(3.5.6) 


Observe  that,  despite  the  restrictions  imposed  on  the  coefficients  of  thermoelastic 
coupling  by  (3.5.1),  the  corrections  to  the  out-of-plane  displacement  and  absolute 
temperature  fields  reamain  coupled  through  (3. 5. 3)2  and  (3.5.6). 

The  initial  conditions  satisfied  by  w  and  s  are,  from  (3.2.8),  (3.2.9)  and 
(3.2.11), 


(3.5.7) 


satisfied  by 

w 

and 

s  are, 

*»(•» 

0 

+, 

II 

V 

on 

2R2, 

w(-, 

0 

+ 

II 

w 

on 

m2, 

T(; 

0 

i 

11 

on 

R\ 

s( 

11 

^ — v 

+ 

0 

h 

on 

R. 

Finally,  from  (3.2.6)  and  (3.2.7),  it  is  assumed  that,  for  each  t  in  2R+,  the  following 
far  field  decay  conditions  hold 


lim  w,i  (x\.  •,?)  =  0  on  1R, 

Xl  — »±oo 


lim  tv#  (-,Z2,<)  =  0  on  2R, 

xj— *±oo 


lim  T(xi,X2,t)  =  0. 

c*+x  *— »±oo 


(3.5.8) 
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In  the  inertia-free  case,  (3.5.2)i  2  are  replaced  by 

a2to,  11  +b]w,22  =  0  on  $/, 

(3.5.9) 

a*W,u  +&rtt,»22  =0  On  ttr. 

Furthermore,  the  jump  condition  (3.5. 3)i  is,  by  virtue  of  (3.4.15),  replaced  by 

[a2tw,i  ]]  =  0  on  /,  (3.5.10) 

and  (3.5. 3)2  is,  from  (3.4.16)  and  (3.5.1),  supplanted  by 

[kT,i  J  +  pvo [c]T  =  pv0( 7r  -  7i)((a2it;,i ))  +  £0s  on  I.  (3.5.11) 

while  (3.5.3)3,4  continue  to  hold.  Finally,  the  linearized  kinetic  relation  (3.5.6)  is, 
upon  referring  to  (3.4.17)  and  (3.5.1),  superceded  by 

i  =  2^«a2u.,.  »+(<>»  + °n  /.  (3.5.12) 

In  the  absence  of  inertial  effects  initial  conditions  cannot  be  given  for  the  incre¬ 
ments  to  the  out-of-plane  displacement  and  velocity  fields  w  and  w\  the  initial 
condition  (3.5.7)33  pertaining  to  T  and  s  still,  however,  continue  to  be  applicable. 
The  decay  conditions  (3.5.8)  also  still  hold. 

3.6.  Normal  mode  analysis  for  a  base  process  involving  a  static 
interface  in  the  absence  of  inertia.  Suppose  that  t>o  in  (3.1.1)  is  zero.  Then 
the  base  process  described  by  (3.1.1)— (3.1.3)  is  a  piecewise  homogeneous  isother¬ 
mal  two-phase  state  involving  a  static  planar  interface.  Recall,  from  Section  2.5, 
that  when  uo  =  0  and  the  kinetic  response  function  V  or  (p  is  continuously  dif¬ 
ferentiable  on  its  domain  of  definition  then  v*  >  0.  Since  u»  =  0  is  ruled  out  by 
whichever  of  (3.1.13)  is  appropriate  and  the  corresponding  expression  (3.4.10)i 
or  (3.4.11)i,  it  is  clear  that — in  the  present  context — u,  >  0.  Consider,  now,  the 


-39- 


initial  boundary  value  problem  composed  by  (3.5.9),  (3.5.2)3)4,  (3.5.10),  (3.5.11), 
(3.5.3)3t4,  (3.5.12),  (3.5.7)3)4  and  (3.5.8).  Note,  since  i»o  =  0,  that  (3.5.11)  and 
(3.5.12)  reduce  with  the  aid  of  (3.1.12)  to 


l*T,i  J  =  pdQ(rji  -  Tfr )s  on  7, 


s  = 


7<  “  7r  „  o...  «  ,  *?!-*? r 


((a  «M»  + 


T  on  7. 


(3.6.1) 


Observe  that,  by  virtue  of  the  linearization,  the  relevant  partial  differential  equa¬ 
tions,  jump  conditions  and  kinetic  relation  are  all  linear  with  constant  coefficients; 

O  i  O 

note,  also,  that  the  domains  17 and  17 q  are  rectangular.  It  is  therefore  possible 
to  find  a  solution  to  the  linearized  partial  differential  equations,  jump  conditions 
and  kinetic  relation  in  the  form 


w(xi,x2,t) 


T(xi,x2it) 


f  Wie+(lX3eiKX3ept  V(xi,x2 ,t)  G  17q  x  IR+, 
1  Wre-^'e******  V(xi,x2 ,<)  6  17J  x  2R+, 

j  0,e+<’lXleiKX3ept  V(zux2,t)  6  17'  *  1R+, 
1  0re~<‘rXleiKX3ept  V(xi,x2,t)  ei0rx  1R+, 


s(x2,t)  =  SeiKX3ept  V(x2,<)  E  Mx  1R+, 


(3.6.2) 


where  the  amplitudes  Wi,  Wr,  0i,  0r  and  5,  wave-numbers  £j,  £r,  C/,  Cr  and  «,  and 
growth-rate  p  are  all  constants.  To  comply  with  the  decay  conditions  (3.5.8)i(3  it 
is  clear  that  &(£/),  &(£r),  &(0)  and  9fc(£r)  must  all  be  positive.  Although  (3.6.2) 
is  not,  in  general,  consistent  with  neither  the  initial  conditions  (3.5.7)3)4  which 
hold  in  the  absence  of  inertial  effects  nor  the  decay  conditions  (3. 5. 8)2,3,  since  <f> 
and  h  are  stipulated  to  be  square  integrable  on  1R,  and  hence  can  be  represented 
as  Fourier  integrals — 


+  OO 

4>{x\,x2)=  J  <j>(x\,K.)etKX3  d.K  V(x1?x2)  €  B2, 

—  OO 


+  OO 

h(x 2)  =  ^  J  h(K)eiKX3  dK 

—  OO 


Vx2  €  IR, 


(3.6.3) 
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it  is  reasonable  to  expect  that  stability  results  can  be  obtained  by  a  normal-mode 
analysis ;  such  an  analysis  entails  substitution  of  (3.6.2)  into  (3.5.9),  (3. 5.2)3^, 
(3.5.10),  (3.5.11),  (3.5.3)3,4  and  (3.5.12)  to  determine  the  growth-rate  p  as  a 
function  of  the  complex  wave-numbers  £i,  £r,  (i,  Cr  and  the  real  wave-number  k. 
If  there  exists  a  complex  growth-rate  p  with  positive  real  part  which  arises  as  a 
solution  to  the  aforementioned  problem  then  the  base  process  will  be  referred  to 
as  linearly  unstable.  Otherwise,  the  base  process  will  be  called  linearly  stable. 

Observe  that  the  amplitudes  0/,  0r  and  S  and  wave-numbers  C/>  Cr  and 
k  must  be  viewed  as  given  for  the  normal  mode  analysis  to  prove  effective  in 
determining  necessary  and  sufficient  conditions,  via  the  analysis  of  a  dispersion 
relation  like  that  performed  in  [9],  for  the  linear  instability  of  the  base  process  with 
respect  to  arbitrary  disturbances  contained  in  the  class  of  perturbations  described 
in  Section  3.2.  Before  proceeding,  note,  from  (3.5.3)4,  that  0/  =  0r  =:  0.  Now, 
substitution  of  (3.6.2)  into  (3.5.9),  (3.5.2)3(4,  (3.5.10),  (3.5.11),  (3.5.3)3,4  and 
(3.5.12)  yields  the  following  relations 


W,  = 
6  = 


u2\k 


aibi(ji  -  7r) 

if  I  |  ,  _  * 

’iMi  fr  — 

Or 


s,  wr  = 


^2M  5  0_  pGq(k,p)s 

drKilr  -  7»)  Vl  -  Vr 


—  |k|,  Cr  =  — 1*|,  0  =  \A2  +p/<*h  Cr  =  V*2  +p/°r,  (3.6.4) 

a\ 


p+  — (i/2|k|  +p0o(«,p))  =  0, 


where  Gq  :  1R  x  (U  —*  Wis  given  by 


G0(k,p) 


_ Mnt  -  nr)2 _ 

ciy/a'fK2  +  atp  +  cry/a2K2  +  arp 


V(k,p)  6  2R  x  (T, 


(3.6.5) 


and  the  constant  v 2  is  defined  as  follows: 


_  aibiarbrjn  -  7r)2 

aibi  +  arbr 


(3.6.6) 


It  is  clear  from  (3.6.4)3i6)7  that,  for  (3.6.2)  to  represent  a  solution  to  (3.5.9), 
(3.5.2)3,4,  (3.5.10),  (3.5.11),  (3.5.3)3,4  and  (3.5.12),  the  amplitude  0  and  the 
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wave-numbers  Q  and  £r  cannot  be  chosen  independently  of  S  and  k.  Hence,  the 
normal  mode  analysis  is  only  of  use  in  analyzing  the  linear  stability  of  the  base 
process  at  hand  with  respect  to  a  certain  class  of  perturbations;  that  is,  it  is  only 
possible — via  this  analysis — to  determine  conditions  necessary  and  sufficient  for 
the  instability  of  the  base  process  with  respect  to  a  proper  subset  of  the  class 
of  perturbations  introduced  in  Section  3.2.  To  achieve  such  results  it  suffices  to 
analyze  the  zero  structure  of  the  dispersion  relation  (3.6.4)8  as  a  function  of  the 
growth-rate  p  for  fixed  values  of  the  wave-number  «  and  the  parameters  7/,  7r, 
t»o,  0o,  <*i,  «r,  61,  K,  <*!,  <*r5  ci,  cr,  p  and  u„.  This  is  done  below. 

To  comply  with  the  restriction  that  and  3 ?(Cr)  are  both  positive,  the 

square  roots  which  appear  in  the  definition  of  Go  are  defined  so  that  for  p  in  1R, 


<x]k2  +  aip  >  0  \J  ot2  k2  +  a  ip  >  0 

a2/c2  +  <*rp  >  0  =»  y/a2K2  +  arp  >  0 

from  which  it  is  clear  that  for  p  in  <F, 


V/c  6  2R , 
V/c  6  1R, 


3?(<*2/c2  +  aip)  >  0 
3i(a2/c2  +  arp)  >  0 


(^Jafn2  +  o/p^ 
31  (^\Jot2K2  +  arp) 


>0 

>0 


V/c  €  2R, 
V/ce  iR. 


(3.6.7) 


(3.6.8) 


Furthermore,  it  is  evident  from  (3.6.8)  that 


£(p)>  0  <=>  3J(G0(/c,p))  >  0  V/cei7.  (3.6.9) 


This  result  shows  that  there  cannot  exist  a  root  p  in  <F  with  3 l(p)  >  0  to  (3.6.4)8 
unless  v*  <  0.  Since  v,  >  0  it  is  clear  that,  at  present,  there  cannot  exist  a  p  in  (T 
to  (3.6.4)8  with  9fc(p)  >  0  for  any  k  in  1R  \  {0}.  Hence,  when  vq  =  0  and  inertial 
effects  are  disregarded  the  base  process  described  in  Section  3.1  is  linearly  stable 
with  respect  to  all  perturbations  within  the  narrowed  set  under  consideration. 

If,  in  place  of  the  foregoing  normal  mode  analysis,  a  full-fledged  Fourier- 
Laplace  transform  analysis  of  (3.5.9),  (3.5.2)3i4,  (3.5.10),  (3.5.11),  (3.5.3)3)4  and 
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(3.5.12)  is  performed,  then  the  narrowing  of  the  class  of  initial  data  necessitated  by 
the  normal  mode  analysis  does  not  occur.  Furthermore,  in  this  case  it  transpires 
that  the  Fourier- Laplace  transform  of  s  can  be  expressed  in  the  form 

h(K)  +  *^H(k,p) 

S(k,P)  =  ’  V(*,p)  6  JZ  x  <F,  (3.6.10) 

where  h  is  the  Fourier  transform  of  h  and,  for  each  (k ,p)  in  JR  x  <F,  H(n,p)  is 
a  functional  of  the  initial  data  77,  w  and  <f>.  From  the  foregoing  discussion  it 
is  apparent  that,  since  v,  >  0  at  present,  there  exist  no  unstable  zeros  of  the 
denominator  of  the  expression  on  the  right-hand-side  of  (3.6.10).  Hence,  when 
the  base  process  is  one  wherein  the  associated  phase  boundary  is  static  prior  to 
the  instant  at  which  the  perturbation  is  imposed  and  inertial  effects  are  ignored, 
it  is  linearly  stable  with  respect  to  all  perturbations  within  the  class  introduced 
in  Section  3.2. 

3.7.  Energy  analysis  for  a  base  process  involving  a  static  interface 
with  inertial  effects  present.  Suppose,  as  in  Section  3.6,  that  Vo  in  (3.1.1) 
is  zero;  the  parameter  v„  is,  as  such,  positive.  Consider,  now,  the  inertial  initial 
boundary  value  problem  formed  by  (3.5.2),  (3.5.3)  and  (3.5.6)-(3.5.8).  Observe 
that,  since  uo  =  0,  (3.5. 3)2  and  (3.5.6)  are  replaced  by  (3.6.1)i  and  (3. 6.1)2,  re¬ 
spectively.  Furthermore,  (3.5. 3)i  simplifies  to  its  inertia-free  counterpart  (3.5.10). 
In  place  of  a  normal  mode  analysis  like  that  performed  in  Section  3.6  an  energy 
analysis  will  be  used  in  this  section  to  show  that,  when  inertial  effects  are  ac¬ 
counted  for  but  Vo  =  0,  the  base  process  described  by  (3.1.1)— (3.1.3)  is  linearly 
stable  with  respect  to  all  perturbations  of  the  type  put  forth  in  Section  3.2.  Pre¬ 
liminary  to  doing  so  define  the  total  energy  S  :  [0,  f,)  by 


£(t)  =  EK(t)  +  Ew{t)  +  ET{t )  Vt  <E  [0 ,<.), 


(3.7.1) 
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where  Ek  :  [0, i«)  — »JR+  is  the  kinetic  energy  given  by 


+oo  +oo 


Ex(t)  =  y  J  J  w(xi,X2ft)dxidx2  Vt€[0,t,),  (3.7.2) 


— oo  — oo 


Ew  ’  [0,  <*)  — is  the  elastic  energy  defined  via 


+oo  0 


Ew(i)  =y J  J {afw,](xi,x2,t)  +  tiw,l(xi,x2,t))dxidx2 


— oo  — OO 
+ oo  -f  oo 


+Y J  J (ar w*(xUX2,t)  +  blw,l(xi,x2,tj)dxidx2 

vt  e 


— oo  0 


(3.7.3) 


and  Et  :  [0,i*)  -*1R+  is  the  thermal  energy  given  by 


+oo  0 


+oo  0 


Er(t)  =y J  J clT2(xi,X2,t)dxidx2  +  ^J  J crT2(xi,x2,t)dx1dx2 


— oo  —OO 
t  +oo  +oo 


— oo  — oo 


+ptj  J  J klT,a(xi,X2,r)T,a  (xi,x2,T)dxidx2dT 


0  — oo  0 
t  +oo  +oo 


+ptj  J  J krT,a(xiyX2,T)T,a(xi,X2,T)dxidx2dT 

Vi  G  [0,  i»).  (3.7.4) 


0  — oo  0 


The  constant  £  which  appears  in  (3.7.2)-(3.7.4)  is  assumed  to  be  positive  and 
carries  units  of  length.  It  is  clear  from  (S.SJ):^^  and  the  stipulated  square 
integrability  of  17, a,  ro  and  <f>  that  £(0)  exists  In  writing  (3.7.2)-(3.7.4)  it  is 
assumed,  however,  that  there  exists  a  positive  time  t,,  which  may  possibly  be  very 
small,  such  that  the  relevant  integrals  exist  on  [0,i,).  A  reasonable  definition  of 
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linear  stability  is,  at  present,  that  £  remain  bounded  on  M+.  A  straightforward 
but  long  calculation  which  makes  use  of  (3.5.8),  (3.6.1)  and,  recalling  the  foregoing 
remarks  regarding  the  coincidence  of  (3.5.3)i  with  (3.5.10)  when  uo  =  0,  show 
that  the  power  £  is  given  by 


+  0O 


—  OO 


(x2,t)dx  2 


Vt  6  (0,t*). 


(3.7.5) 


Since  uo  =  0  at  present,  £(t)  <  0  for  all  t  in  [0,t*);  xmder  these  circumstances  the 
interval  over  which  £  is  defined  can  be  extended  incrementally  to  JR.+  leading  to 
the  following  inequality: 

£(t)  <0  V*  6  R+.  (3.7.6) 

Evidently,  then,  by  the  definition  of  linear  stability  given  above,  the  base  process 
at  hand  is  stable  with  respect  to  all  perturbations  under  consideration  if  the 
associated  phase  boundary  is  static  prior  to  the  instant  at  which  the  perturbation 
is  imposed  and  inertial  effects  are  accounted  for. 

Note  that  a  normal  mode  analysis  akin  to  that  performed  in  Section  3.6 
produces  the  following  dispersion  relation 

P+  —  (Fo(k,p)+pG0(k,p))  =  0,  (3.7.7) 


where  Fq  :  1R  x  (T  — ►  <F  is  defined  by 


Fq(k,p)  -  TdV*  VKp)6Rx( V, 

aibiy/K 2  +p2/tf  +  Or  bry/K2  +  p2 /b2r 

(3.7.8) 

and  Go  is  as  defined  in  (3.6.5).  A  study  of  this  dispersion  relation  allows  the 
recovery  of  the  results  obtained  by  the  foregoing  energy  analysis. 

The  combined  results  of  this  and  the  proceeding  section  are  consistent  with 
those  presented  by  Fried  [9]  in  the  purely  mechanical  analogue  of  the  problem 
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considered  here.  Hence,  when  vo  =  0,  the  inclusion  of  thermal  effects  does  not  al¬ 
ter  the  linear  stability  of  the  base  state  (3.1.1)-(3.1.3)  from  its  obvious  mechanical 
analogue. 

3.8.  Normal  mode  analysis  for  a  base  process  involving  a  moving 
interface  with  or  without  inertial  effects.  Suppose  that  v0  in  (3.1.1)  is 
positive.  Consider  now  both  the  inertia-free  intial  value  problem  consisting  of 
(3.5.9),  (3.5.2)3,4,  (3.5.10),  (3.5.11),  (3.5.3)M,  (3.5.11),  (3.5.7)3>4  and  (3.5.8)  and 
the  inertial  initial  value  problem  comprised  by  (3.5.2),  (3.5.3)  and  (3.5.6)-(3.5.8). 
Note  that,  in  both  of  these  cases,  (3.1.5)3,  (3.1.11)  and  the  assumed  positivity  of 
v0  imply  that  =  0.  Hence,  in  the  inertia-free  case,  (3.5.11)  and  (3.5.12)  simplify 
as  shown  below 


V°Ty\  ]  +  pv0[c]T  =  pv o(7r  -  7/)((a2iv,i ))  on  /, 

g  _  1 1 — lL([a y  +  VdT  on 

v* 

while,  in  the  inertial  case,  (3.5.3)2  and  (3.5.6)  become 

lkT, i  ]  +  pv0[c\T  =  pv0(7r  -  7/)(((fl2  -  t>o)«M  ))  on  /, 
s  =  — — —{((a2-vl)w,x))  +  vtT  on  I. 


Next,  a  normal  mode  analysis  analogous  to  those  undertaken  in  Section  3.6  and 
[9]  will  be  performed  based  upon  the  following  representation  of  a  solution  to  the 
relevant  partial  differential  equations,  jump  conditions  and  kinetic  relation: 


w(xi,x2,t) 


T(xi,x2,t) 


W|C+(*.-voO*ic;«Scp<  V(xi,x2)  €  fi[,  t  €  54, 

V(x1,x2)€/7[,  <  €  54, 

e,e+(G-»o«)*iciK*aeFt  V(x,,x2)  €  fi[,  <  €  54, 

0rC— «r-vo t)xleiKziept  v(x,,x2)e^r,  <6  5 4, 


s(x2,<)  =  Se*KXlept  V(x2,<)  6  R  x  54- 


(3.8.3) 
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Substitution  of  (3.8.3)  into  the  equations  appropriate  to  the  equations  appropriate 
to  the  inertia-free  and  inertial  cases  gives,  respectively,  the  following  dispersion 


relations  2 

p+^-!S(i  +  ^g(/c,p))  =  o, 
f\±p)(1+vtG(Kp))  =  0f 

V , 

where  G  :  1R  x  <T->  (Tand  F  :  1R  x  <F -*• (Tare  given  by 

GU.p)  - - - ; - 7  V(/c,p)  £  x  <r, 

V  C/  -  cr  +  ctgi{K,p)  +  crgr(ti,p) 

Fu,p)  =  to  -  +  ">■!)  V(K,p)  6  iR  X  <r, 

V  '  /l(*,p) +  /.-(*,  P) 

where  gt:Rx  <F->  (Tand  gr  :  2R  x  (T  <F  are  given  by 

g/(K,p)  =  ^1  +  ^(ai  *2  +  QiP)  V(K’P)  €  2R  x  <F, 
5r(«,p)  =  +  ^(«r«2  +  arp)  V(k,p)  eJRxW, 

and  fr.  Rx  (E-*  <T and  fr:Rx  <E -*  W are  defined  via 

/,(«,p)  =  -  vl)b]K2  +  a?P2  V(«,p)  E  Rx(E, 

/r(«,p)  =  ^/(a?  ~  vl)b2rK2  +  a?P2  V(/c,p)  E  RxW. 

When  vrf  =  0  the  dispersion  relations  in  (3.8.4)  reduce  to 

p  +  i^Uo, 
p+^  =  0, 


(3.8.4) 


(3.8.5) 


(3.8.6) 


(3.8.7) 


(3.8.8) 


Observe  that  (3.8.8)i  and  (3.8.8)2  are  structurally  identical  to  the  inertia-free 
and  inertial  dispersion  relations  obtained  by  Fried  [9]  in  the  purely  mechanical 
analogue  of  the  investigation  at  hand.  Hence,  if  the  kinetic  response  function  V 


-47- 


or  <p  is  chosen  so  that  =  0  then  the  linear  stability  of  the  base  process  at  hand 
remains  unaltered  from  that  of  its  purely  mechanical  analogue  by  the  presence  of 
thermal  effects;  specifically,  when  v#  =  0,  the  linear  stability  of  the  base  process 
(3.1.1)-(3.1.3)  is  determined  entirely  by  the  sign  of  v,.  That  is,  vm  <  0  is  a 
necessary  condition  for  the  base  process  to  be  linearly  unstable  with  respect  to 
any  perturbation  of  the  type  introduced  in  Section  3.2  and,  further,  <  0  is  a 
sufficient  condition  for  the  base  process  to  be  linearly  unstable  with  respect  to  all 
but  a  small  class  of  very  special  initial  disturbances  contained  within  the  full  set 
under  consideration.  Note,  in  particular,  that  v#  =  0  if  either  V  depends  only 
on  9  through  the  ratio  f/8  or  Cp  is  independent  of  9.  Assume,  henceforth,  that 
v*  ^  0. 

The  branches  of  the  square  roots  which  define  gi  and  gr  are  chosen  so  that, 
for  p  in  JR, 


1  +  A(a2K2  +  atip)  >  0 
v0 

1  +  4(«r«2  +  <*rP)  >  0 
0 


gi(K,p)  >0  V/c  G  IR, 
<7r(«,p)  >0  V/c  €  1R, 


from  which  it  is  clear  that,  for  p  in  <F, 


(3.8.9) 


&(1  +  4(a?/c2  +  «/p))  >  0 

£(1  +  A(cxIk2  +  arp))  >  0 
0 


&(£/(«>?))  >0  Vic  6  JR, 

(3.8.10) 

^Csv^p))  >0  V/c  g  ir. 


Evidently,  then,  (3.8.10)  and  (3.8.5)i  yield  the  following  result: 


5R(p)  >  0  ■$=>  R(v*G(/c,p))  >0  V/c  G  JR.  (3.8.11) 


The  square  roots  which  appear  in  fi  and  fT  axe  defined  via  the  principal  branch 
of  the  complex  logarithm.  It  is,  therefore,  clear  that 

£(p)>0  4=*  B(F(/c,p))>  0  V/c  G  IR.  (3.8.12) 


-48- 


An  immediate  consequence  of  (3.8.11)  and  (3.8.12)  is  that  one  or  both  of  the 
parameters  vm  or  v$  must  be  negative  in  order  for  a  root  p  in  <F of  either  (3.8.4)i  or 
(3.8.4)2  to  have  a  positive  real  part.  It  is  also  obvious,  from  (3.8.11)  and  (3.8.12), 
that  if  both  vm  <  0  and  v$  <  0  then  there  exists  a  root  p  in  (T  with  Bt(p)  >  0  to 
both  of  (3.8.4)  regardless  of  the  value  of  the  wave  number  k  in  1R.  \  {0}.  A  more 
subtle  condition  sufficient  for  the  existence  of  a  root  p  in  <F  to  either  of  (3.8.4) 
occurs  under  the  assumption  that  u*  >  0,  v#  <  0  and  vm\v$\/ci  >  1.  Specifically, 
when  v,  >  0,  v#  <  0  and  vm\v#\/ci  >  1  it  is,  then,  possible  to  show  that  there 
exists  a  root  p  in  (T  to  both  of  (3.8.4)  provided  the  wave-number  k  in  JR  \  {0}  is 
sufficiently  small  so  that  the  inequality 


Cl-Cr  fl(K,p)  fr(li,p)  2v.\v*\ 
Cl  +  cr  1  +  ^  1  +  ^  Cl  +cr 


(3.8.13) 


holds.  A  similar  condition  which  guarantess  the  existence  of  a  root  p  in  (F  to  either 
of  (3.8.4)  occurs  under  the  assumption  that  v,  <  0,  v#  >  0  and  \v+\v$/ci  <  1. 
In  this  case  there  always  exists  such  a  root  to  either  of  (3.8.4)  as  long  as  the 
wave-number  k  is  sufficiently  large  so  that  the  following  inequality  is  satisfied: 


Cl-Cr  fl(K,p )  fr(K,p)  2\v.\v* 
Cl  +  cr  +  1  +  ^  +  1  +  ^  Cl  +  cr 


(3.8.14) 


The  foregoing  discussion  shows  that,  unlike  the  purely  mechanical  process 
investigated  in  [9],  the  present  context  is  not,  when  Vo  >  0,  amenable  to  the 
statement  of  necessary  and  sufficient  conditions  for  the  linear  instability  of  the 
base  process  at  hand.  The  sufficient  conditions  which  have  been  presented  above 
are,  however,  of  interest. 

3.9  Conclusion.  In  {9]  it  is  demonstrated  that  when  the  purely  mechanical 
analogue  of  the  parameter  v,  is  positive  the  appropriate  purely  mechanical  version 
of  the  base  process  considered  here  is  linearly  stable  with  respect  to  all  perturba¬ 
tions  which  axe  considered  in  that  context.  The  last  of  the  conditions  sufficient 
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for  the  linear  instability  of  the  thermoelastic  base  process  (3.1.1)-(3.1.3),  viz., 

V,V0 


vm  >0,  v#  <  0, 


ci 


>i,  ^  +  +  <  Mss!,  (3.9.i) 


c/  +  cr  1  + 


ci 


i  +  S 


Cl  +  cr 


where  the  parameters  u,,  v$,  c/  and  cr  are  as  defined  in  (3.4.10)i  or  (3.4.11)i, 
(3.4.10)2  or  (3.4.11)2,  (3. 5.5)4  and  (3. 5. 5)3,  respectively,  and  the  functions  //  and 
fr  are  as  defined  in  (3.8.7),  is,  hence,  arguably  the  most  interesting  of  the  three 
which  are  presented.  It  exposes  what  might  be  described  as  a  competition  between 
mechanically  stabilizing  and  thermally  destabilizing  effects  and  an  explicit  depen¬ 
dence  of  growth-rate  upon  wave-number.  Significantly  in  these  circumstances, 
it  is  the  low  wave-numbers  (that  is,  long  waves)  with  respect  to  which  the  base 
process  is  linearly  unstable.  Under  conditions  consistent  with  (3.9.1),  a  moving 
planar  phase  boundary,  therefore,  tends  to  prefer  a  highly  wrinkled — i.e.,  plate¬ 
like  or  dendritic — morphology.  Instability  of  this  variety  is  also  found  in  models 
for  dendritic  crystal  growth  and  solidification  (see  [13-14]  and  [16]). 

In  analogy  to  [9]  where  the  physical  plausibility  of  a  purely  mechanical  kinetic 
response  function  for  which  the  parameter  analogous  to  u*  can  be  negative  is 
addressed,  it  is  now  natural  to  consider  the  question  of  whether  it  is  physically 
reasonable  for  a  kinetic  response  function  to  depend  monotonically  on  its  first 
argument — so  that  vm  is  always  positive  and  the  related  purely  mechanical  process 
is  linearly  stable — but  non-monotonically  on  its  second  argument,  in  which  case 
vj  may  be  negtive.  The  experimental  work  of  Clapp  &  Yu  [5]  which  studies,  in 
part,  the  dependence  upon  temperature  of  transformation  kinetics  in  a  particular 
alloy  capable  of  sustaining  displacive  solid-solid  phase  transformations  indicates 
that  the  role  of  temperature  in  such  kinetics  is  very  complicated.  In  fact,  despite 
what  appears  to  be  a  very  careful  experimental  procedure  and  analysis,  Clapp 
&  Yu  [5]  observe  a  severe  scatter  in  the  data  which  measure  the  dependence  of 
phase  boundary  velocity  upon  temperature.  This  scatter  indicates  there  may  not 
be  a  simple  functional  dependence  of  interface  normal  velocity  upon  temperature. 
With  regard  to  the  issue  at  hand,  these  experimental  results  seem  to  indicate  that, 
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if  a  kinetic  relation  of  the  form  (2.5.1)  or  (2.5.2)  is  insisted  upon,  monotonicity 
of  a  kinetic  response  function  V($,  •)  for  fixed  #  in  2R  or  </?(V,  •)  for  fixed  V  in  1R 
may  be  the  exception  rather  than  the  rule. 
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r(  7i^o) 


Figure  1:  Graph  of  the  shear  stress  response  function  r(-,  60)  for 
fixed  $o  in  ( Om,0M )• 
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function  V(-,0o)  for  fixed  0o  in  (0m,0M). 
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